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I.  INTRODUCTION 


Nonuniform  metallic  waveguide  transitions  have  been  the  subject  of 
many  theoretical  investigations.  A  general  solution  was  first  given  by 
Stevenson,*  who  expanded  the  field  intensities  into  a  series  of  cross- 
sectional  wave  functions.  Later  using  the  same  approach,  Schelkunof f , 2 
Reiter3  and  Katzenelenbaum14  independently  derived  the  generalized  tele¬ 
graphist's  equations;  thus  describing  waveguide  transitions  as  an 
infinite  set  of  coupled  transmission  lines.  Solymar^  transformed  these 
into  a  set  of  differential  equations  for  the  amplitudes  of  the  forward 
and  backward  traveling  waves.  The  power  in  the  equations  these  men 
derived  remained  dormant  until  high  speed  computers  could  be  implemented 
to  solve  them. 

The  need  to  efficiently  design  nonuniform  waveguide  transitions 
exists  in  almost  every  facet  of  electromagnetics  engineering.  These 
devices  are  designed  to  maximize  power  transfer  from  one  size  waveguide 
to  another  and  are  called  impedance  matches  or  transformers.  Generally, 
a  costly,  time-consuming  "build  and  test”  method  is  used  to  optimize  the 
power  reflection  and  transmission  characteristics  of  impedance  matches. 
Numerical  methods  can  be  used  to  replace  this  procedure  by  a  practical 
computerized  design  tool. 

The  purpose  of  this  work  is  to  outline  the  theoretical  and  numer¬ 
ical  aspects  of  developing  a  nonuniform  waveguide  transition  design 
tool,  and  use  them  to  write  a  computer  program  that  models  dominant  mode 
rectangular  and  ridged  waveguide  tapers. 


Theoretical  and  numerical  aspects  of  design  tool  development  are 
discussed  in  Sections  II  and  III,  respectively.  Section  II  summarizes 
the  work  of  Reiter  and  Solymar  and  outlines  a  method  by  which  their 
formulae  can  be  used  to  obtain  a  transition's  scattering  matrix.  In 
Section  III,  aspects  of  numerically  modeling  the  theoretical  formulae  of 
Section  II  are  presented  in  parallel  with  the  development  of  the  domi¬ 
nant  mode  computer  program  for  double-ridged  waveguide  transitions.  The 
final  section  shows  that  computed  standing  wave  ratios  (VSWR)  of  ridged 
and  unridged  transitions  agree  well  with  experimentally  measured  values. 


II.  A  THEORETICAL  ANALYSIS  OF  THE  NONUNIFORM  WAVEGUIDE  TRANSITION 


Presented  herein  is  a  formulation  that  shows  how  the  scattering 
matrix  of  a  tapered  waveguide  can  be  obtained  from  a  set  of  transmission 
line  equations  that  model  it.  A  review  of  the  work  performed  by  Reiter3 
bridges  the  gap  between  the  Maxwell  equations  and  the  infinite  set  of 
voltage-current  differential  equations  that  model  a  waveguide.  A  sum¬ 
mary  of  Solymar's^  work  shows  how  these  differential  equations  were 
modified  in  order  to  describe  the  transition  in  terms  of  forward  and 
backward  traveling  waves.  The  transmission  matrix  obtained  by  solving 
these  traveling-wave  equations  is  algebraically  transformed  into  the 
scattering  matrix  of  the  taper.  The  formulation  begins  with  a  concise 
mathematical  and  physical  description  of  the  taper. 

A.  A  Description  of  the  Problem 

The  problem  is  to  determine  the  scattering  matrix  of  the  tapered 
waveguide  transition  depicted  in  Fig.  1.  The  transition  consists  of  a 
tube  bounded  by  a  conducting  surface  such  that  a  plane  perpendicular  to 
the  z-axis  cuts  this  surface  in  a  closed  curve  C(x,y,z).  The  region  cut 
out  by  C  is  the  cross  section  of  waveguide  denoted  by  A.  Both  A  and  C 

A  A 

are  considered  to  be  continuous  functions  of  z.  The  symbols  n  and  z 
denote  the  normal  unit  vectors  to  C  and  A,  respectively.  s  is  along  C 

A  A 

and  is  perpendicular  to  both  n  and  z.  0(x,y,z)  is  defined  by  the  angle 

A  A  A 

between  z  and  a  line  in  the  n-z  plane  which  is  tangent  to  the  waveguide 
wall  at  A.  Waveguides  Guide  I  and  Guide  2  are  fed  by  some  linear  combi¬ 
nation  of  pure  modes  and  are  assumed  to  extend  infinitely  beyond  z  *  0 


Fig.  1.  Physical  and  mathematical  description 
of  a  nonuniforra  waveguide  transition. 

and  z  *  L.  These  definitions  are  next  used  to  model  a  transition  by 
current-voltage  differential  equations. 

B .  The  Generalized  Telegraphist’s  Equations 

This  review  of  Reiter's  work  shows  that  an  electromagnetic  wave 
propagating  within  the  confines  of  a  waveguide  is  equivalent  to  a  system 
of  coupled  transmission  lines.  Reiter  derived  this  system  of  transmis¬ 
sion  line  equations  by  enforcing  waveguide  wall  boundary  conditions  upon 


z)  'ep 


7t  is  the  transverse  gradient  operator  and  ip  is  the  mode’s  longitudinal 

field  function.  <p  is  proportional  to  H_  and  E  for  TE  and  TM  modes, 

P  z  z 

respectively.  Et  and  Ht  are  also  functions  of  ip. 

The  waveguide  wall  boundary  conditions  are  enforced  upon  the 
Maxwell  equations  through  ip^  and  h^.  ip^  and  hp  are  obtained  by  solving 
the  transverse  Helmholtz  wave  equation. 


(V  2  +  h2)  ip  -  0 
v  t  p '  p 


The  boundary  conditions  on  ip  are  ip  =0  for  TM  modes  and  9ip  / 9n  =  0  for 

P  P  P 

TE  modes.  There  are  an  infinite  number  of  field  configurations  ip  which 
satisfy  each  of  these  boundary  conditions.  Each  of  these  field  config¬ 


urations  has  a  corresponding  eigenvalue  h.  As  stated  in  Eqs.  2  and  3,  a 
linear  combination  of  these  modes  is  sufficient  to  describe  the  electric 
and  magnetic  fields  of  a  bounded  wave. 


By  using  various  mathematical  identities  in  conjunction  with  Eqs. 
1  through  6,  the  Maxwell  curl  equations, 


+  •+■  -*• 
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7  x  h  «=  jwe  E 


$  •  E  =  0 


7  •  H  *  0 
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are  transformed  into  an  equivalent  system  of  coupled  transmission  line 
equations  given  by 


The  subscripts  i  and  p  denote  arbitrary  modes.  and  ic^  are  the  mode 
propagation  constant  and  wave  impedance,  respectively.  The  Tp^  and  Tj^ 
voltage  and  current  transfer  coefficients  are  listed  in  Appendix  A. 
Equation  10  is  the  lossless  form  of  the  Generalized  Telegraphist's 
Equation.  It  shows  that  a  waveguide  is  equivalent  to  a  coupled  system 
of  transmission  lines.  This  system  of  voltage-current  equations  can  be 
rewritten  in  a  form  which  describes  propagation  within  a  taper  in  terms 
of  forward  and  backward  waves. 


C.  The  Normal  Mode  Equations 

By  choosing  an  appropriate  linear  combination  of  the  transmission 
line  voltages  and  currents,  Eq.  10  can  be  transformed  into  a  system  of 
differential  equations  written  in  terms  of  the  amplitudes  of  forward  and 
backward  traveling  waves.  It  is  customary  to  choose  a  linear  combina¬ 
tion  that  makes  the  magnitude  squared  of  the  mode  amplitudes  propor¬ 
tional  to  the  time  average  power  carried. 

The  linear  combination  of  mode  voltages  and  currents  given  by 
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is  defined  such  that  the  net  time  average  power  flowing  in  the  +z  direc¬ 


tion  at  any  point  z  is  given  by 


1  ,+  ■*■* 
j  Re  f  E  x  H  da 

2  A 


1  c  * 

~  Re  l  V  I 

2  P  P 
P 


l\ |a+ 

1  P 
D  L 


(  z  )  I  2  -  |  Ap  (  Z  )  |  2 


Equation  11  can  be  inverted  to  give  V  and  I  as 


Vp  -  (2  "p)1/2  (*p  +  *p) 


I  .  (2/«  )UZ  (A+  -  A") 
p  ^  ^P  P 


When  Eq.  13  is  substituted  into  Eq .  10,  Solymar’s  normal  mode  form  of 


the  Generalized  Telegraphist's  Equations  results: 
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and  are  the  forward  and  backward  coupling  coefficients  given  in 
Appendix  A.  A^  and  A^  are  the  amplitudes  of  the  forward  and  backward 
traveling  waves,  respectively. 

These  normal  mode  equations  reveal  much  about  how  waveguide  modes 

propagate  and  interact.  8^  is  the  mode  propagation  constant.  An  axial 

change  in  the  waveguide  impedance  ic  causes  mode  reflection.  The  S* 

i  ip 

coefficients  are  responsible  for  self  and  intermode  coupling.  They 
depend  on  the  boundary  fields  and  cutoff  frequencies  of  the  i*-*1  and  pc^ 
modes  and  may  be  interpreted  as  arising  directly  from  geometric  effects. 

As  seen  in  Appendix  A,  it  is  convenient  to  employ  the  convention 
of  enclosing  TM  and  TE  mode  subscripts  with  parentheses  and  brackets, 
respectively  (e.g.,  TMqj^  and  TEfQi])*  When  this  notation  is  applied 
to  the  Helmholtz  Wave  Equation,  the  coupling  coefficient  and  mode  field 
solutions  are  written  for  TM  modes  as 
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To  summarize,  Maxwell's  curl  equations  have  been  transformed  into 
normal  mode  equations.  Reiter's  work  made  it  possible  to  represent  a 
waveguide  by  an  equivalent  system  of  coupled  transmission  lines.  An 
appropriate  linear  combination  of  transmission  line  voltages  and  cur¬ 
rents  resulted  in  a  system  of  equations  that  describe  energy  propagation 
in  waveguides  in  terms  of  forward  and  backward  traveling  waves.  These 
normal  mode  equations  are  next  used  to  obtain  the  scattering  matrix  of  a 
tapered  waveguide. 
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D.  A  Scattering  Matrix  Formulation 


The  method  used  to  compute  a  taper's  scattering  matrix  is  pre¬ 


sented  in  three  parts:  1)  convert  the  normal  mode  equations  from  com¬ 


plex  to  real  form,  2)  formulate  a  transmission  matrix  T  and  3)  algebrai¬ 


cally  transform  T  into  the  scattering  matrix  S.  The  first  step  in  this 


formulation  can  be  bypassed  if  a  complex  variable  differential  equation 


solving  routine  is  available. 


1 .  Converting  Complex  Normal  Mode  Equations  into  Real  Ones 


Since  the  differential  equation  solver  used  in  this  work  deals 


with  real  variables,  it  is  necessary  to  convert  the  complex  normal  mode 


equations  into  real  ones.  Equation  14  can  be  written  as 


Both  A  and  M  can  be  separated  into  real  and  imaginary  parts, 


*d t  *±r  x  .t±i 
A  *  A  +  jA 


=  =r  =i 
H  ■  M  +  1M 
pq  pq  ■  pq 


where  the  r  and  i  superscripts  denote  real  and  imaginary  parts,  respec¬ 


tively.  Substituting  Eq .  18  into  Eq.  17,  carrying  out  the  multiplica¬ 


tion  and  separating  real  and  imaginary  parts  yields 
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(19) 


In  order  to  distinguish  the  complex  normal  mode  equations  from  the  real 
ones,  Eqs.  17  and  19  are  rewritten  as 


dz 


x  = 


Mx 


and 


d  -*■ 


d7  y 


-  By 


(20) 


respectively.  For  N  modes,  the  length  of  x  (complex)  is  2N  and  the 

a  = 

length  of  y  (real)  is  4N.  The  dimensions  of  M  and  B  are  2N  x  2N  and 
4N  x  AN . 

Equation  20  is  a  real  matrix  form  of  the  normal  mode  differential 
equations  used  to  model  a  waveguide  taper.  The  2N  coupled  complex 
differential  equations  (2  equations  per  mode,  one  for  each  of  the  for¬ 
ward  and  backward  waves)  have  been  changed  into  4N  real  equations.  The 
mathematical  model  has  been  reduced  to  a  form  suitable  for  numerical 
solution.  The  next  step  is  to  obtain  solutions  and  transform  them  into 
a  transmission  matrix  for  the  tapered  waveguide  section. 
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2 .  Formulating  a  Transmission  Matrix 

The  transmission  matrix  is  obtained  in  the  following  manner:  1) 
use  orthogonal  initial  condition  vectors  to  solve  the  normal  mode  matrix 
differential  equation  and  2)  algebraically  convert  initial  condition  and 

s 

solution  vectors  into  T.  Figure  2  shows  a  two-port  device  for  which  the 
incident  and  transmitted  waves  are  expressed  in  terms  of  normal  mode 
amplitudes . 


Fig.  2.  A  signal  flow  diagram  expressing  incident 
and  transmitted  waves  in  terms  of  normal 
mode  amplitudes  at  z  =  0  and  z  =  L. 

For  this  two-port  device,  the  normal  mode  differential  equations 
depict  a  two-point  boundary  value  problem  (BVP).  The  known  (incident) 
and  unknown  (transmitted)  signals  exist  at  both  ports.  Standard  differ¬ 
ential  equation  solving  routines  solve  initial  value  problems  (unknowns 
at  one  boundary,  knowns  at  the  other).  By  carefully  choosing  the  ini¬ 
tial  condition  vectors,  these  routines  can  be  used  to  solve  the  BVP. 


^  + 


a.  An  Orthogonal  Set  of  Initial  Condition  Vectors 

Any  set  of  mode  amplitude  initial  condition  vectors  can  be  used  to 
obtain  the  transmission  matrix;  however,  it  is  mathematically  convenient 
to  choose  the  orthogonal  set  given  by 
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As  can  be  seen  from  Eq .  19,  the  first  N  initial  condition  vectors  corre¬ 
spond  to  incident  modes  (at  z  =  0)  having  unit  amplitude  and  zero 
phase.  Equation  21  shows  that  for  these  modes,  the  ith  element  of  y  (0) 
is  one.  Likewise,  the  second  N  vectors  correspond  to  reflected  waves 
(at  z  =  0)  with  unit  amplitude  and  zero  phase.  By  using  these  initial 


< 
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conditions  to  solve  Eq .  20,  one  obtains  the  following  set  of  linearly 


independent  solution  vectors  at  z  =  L, 


y^L),  y2(L) 


y2N-l(L^’  y2N^L) 


As  a  matter  of  clarification,  there  are  N  modes,  each  having 


forward  and  backward  waves  with  real  and  imaginary  parts.  This  makes 


the  length  of  the  vectors  in  Eqs.  21  and  22  equal  to  4N.  Since  each 


mode  has  an  initial  condition  on  its  forward  and  backward  components, 


there  are  2N  initial  condition  and  solution  vectors.  y(0)  and  y(L)  are 


next  combined  to  obtain  T. 


b.  The  Transmission  Matrix 


The  transmission  matrix  is  constructed  by  algebraiclly  joining 


linear  combinations  of  the  initial  condition  and  solution  vectors.  This 


process  begins  by  transforming  these  vectors  back  into  their  complex 


form;  hence,  there  are  2N  initial  condition  and  solution  vectors  each 


containing  2N  elements.  To  denote  this  change,  the  notation  of  Eas.  21 


and  22  is  changed  to 


x  2  (0)  ,  x2(0)  ...»  x2N-1^’  *2N^ 


The  first  N  solution  vectors  correspond  to  the  transmitted  portion  of 
the  forward  waves.  The  second  N  solution  vectors  represent  the 
backward  waves  that  would  be  needed  at  2  =  L  to  realize  the  ini¬ 
tial  conditions  on  the  backward  waves  at  z  =  0. 
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x^(L),  x^(L)  •  •  •  )  X2N-1  •  x2n  ^ 


respectively.  A  general  solution  and  initial  condition  vector  can  be 
written  as  a  linear  combination  of  the  vectors  in  Eqs.  23  and  24, 


x(0)  =  l  C  x  (0) 
p=l  P  P 


x(L)  =  l  C  x  (L) 
p=l  P  P 


Equations  25  and  26  may  be  written  in  matrix  notation  as 


where  the  columns  of  the  complex  matrices  T  and  U  are  made  up  of  the 
solution  and  initial  condition  vectors,  respectively.  C  is  a  vector 
made  up  of  the  Cp  coefficients.  Changing  y(0)  of  Eq .  21  into  its  com'- 
plex  form  x(0)  of  Eq.  27  shows  that  U  is  the  identity  matrix 


U  =  I 


Solving  Eq .  27  for  C  gives 
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U  x(0) 

—  1  > 

I  x(0) 

x(0) 


(30) 


Substituting  C  into  Eq.  28  yields 


x(L) 


T  x(0) 


(31) 


T  relates  the  mode  amplitudes  at  z  ■  L  to  those  at  z  ■  0  and  is  identi¬ 
cally  the  transmission  matrix.  It  has  been  constructed  by  using  matrix 
algebra  to  properly  combine  an  orthogonal  set  of  initial  condition 
vectors  for  the  taper’s  normal  mode  equations.  The  scattering  matrix 
can  now  be  determined  by  viewing  Eq.  31  in  terms  of  incident  and 
reflected  waves. 


3.  The  Scattering  Matrix 


The  scattering  matrix  is  algebraically  obtained  from  the  transmis¬ 
sion  matrix  by  changing  from  the  notation  of  forward  and  backward  trav¬ 
eling  waves  to  that  of  incident  and  reflected  signals.  Figure  3  shows 
forward  mode  amplitudes  at  z  **  0  and  z  -  L  with  the  labels  al  and  b2, 
respectively.  Likewise,  the  backward  modes  at  these  planes  are  labeled 
$1  and  a2.  The  subscripts  1  through  N  represent  the  N  propagating 
modes.  In  this  signal  flow  notation,  "a"  and  "b"  represent  the  incident 
and  reflected  components  of  a  mode's  energy.  For  example,  when  the  TEjq 
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Fig.  3.  Scattering  signal  flow  diagram  for 
a  nonuniform  waveguide  transition. 

mode  is  fed  into  a  transition  from  both  ends;  al^  and  a2^  are  propor¬ 
tional  to  the  incident  power  at  ports  1  and  2.  bl^  and  b2j  are  likewise 
proportional  to  the  reflected  power. 

x  E 

Obtaining  S  from  T  is  a  matter  of  algebraically  transforming  Eq. 
31  with  the  notation  of  Fig.  3  incorporated  into  it.  Replacing  x  in  Eq. 
31  by  a  and  b  yields 

(32) 

The  scattering  matrix  relates  incident  signals  to  reflected  ones.  It 
can  be  obtained  by  solving  Eq .  32  for  bl  and  b2 , 


Writing  this  in  matrix  form 


and  reducing  it  to  a  single  equation  yields 


(34) 


(35) 


where  S  is  identically  the  scattering  matrix  of  the  tapered  transition. 

In  summary.  Maxwell's  equations  as  they  apply  to  a  nonuniform 
waveguide  transition,  have  been  solved  to  obtain  a  multimode  scattering 
matrix.  Reiter's  transmission  line  model  of  a  waveguide  was  used  as  the 
starting  point  for  developing  a  set  of  normal  mode  equations.  These 
were  obtained  by  writing  linear  combinations  of  the  transmission  line 
voltages  and  currents  that  defined  the  amplitudes  of  forward  and  back¬ 
ward  traveling  waves.  The  transmission  matrix  was  expressed  as  an 
algebraic  combination  of  the  initial  condition  and  solution  vectors  of 
these  traveling  wave  equations.  Finally,  the  signal  flow  notation  of 
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incident  and  reflected  waves  was  used  to  transform  T  into  S.  The  next 
section  shows  that  with  current  numerical  methods,  this  formulation  can 
be  used  to  obtain  the  scattering  matrix  of  transitions  in  rectangular 
and  double-ridged  waveguides. 
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III.  NUMERICAL  DESIGN  TOOL  DEVELOPMENT  FOR  DOUBLE  RIDGED  WAVEGUIDES 


Two  major  aspects  of  writing  a  computer  program  that  is  capable  of 
modeling  an  arbitrarily  shaped  waveguide  transition  are:  1)  ascertain 
the  axial  dependency  of  the  eigenvalues  and  coupling  coefficients  for 
each  mode  and  2)  obtain  the  transition  scattering  matrix  by  solving  the 
coupled  system  of  differential  equations.  To  illustrate  the  practical¬ 
ity  of  implementing  the  technique  presented  in  Section  II,  a  code  was 
developed  (Appendix  B)  that  models  continuous  symmetrical  double-ridged 
waveguide  tapers  operating  in  the  TE^q  mode.  This  section  gives  a 
detailed  description  of  how  the  finite  difference  method  was  used  to 
compute  coupling  coefficients.  It  also  shows  how  the  single  mode 
coupled  differential  equation  solutions  are  transformed  into  S. 

A .  Numerically  Obtaining  the  Coupling  Coefficients 

Computing  the  axial  dependency  of  the  coupling  coefficients  is 
described  in  three  parts:  1)  numerically  solving  the  transverse 
Helmholtz  Wave  Equation,  2)  showing  how  these  solutions  are  used  to 
obtain  the  coefficients  and  3)  using  a  cubic  spline  to  approximate  an 
axially  discrete  coupling  coefficient  profile  by  a  continuous  one.  For 
some  waveguide  cross  sections,  solutions  to  the  Helmholtz  Wave  Equation 
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t  p  pTp 


(36) 


can 


be  expressed  in  analytical  form.  In  general,  however,  the  p 


th 


Sylvester's6  classic  finite  difference  scheme  was  used  to  find  the 
TEjq  mode  eigenvalue  and  eigenfunction  of  a  double— ridged  waveguide.  A 
discussion  of  methods  that  can  be  used  to  analyze  other  geometries  is 
given  by  Davies7  and  Ng.8  The  following  is  a  brief  summary  of  the  way 
Sylvester's  method  describes  the  cross-section  shape,  the  Helmholtz 
equation,  its  boundary  conditions  and  a  solution  procedure  to  the  com¬ 
puter. 

The  finite-difference  method  uses  a  rectangular  mesh  to  mathemat¬ 
ically  model  the  geometry  of  Fig.  4.  In  order  to  keep  the  analysis 
simple,  a  square  mesh  was  used  and  symmetries  of  the  TEjq  mode  were 
utilized.  The  dotted  vertical  and  horizontal  lines  represent  planes  of 
even  and  odd  H.,  symmetry,  where  the  z-axis  is  into  the  page.  Figure  5 
shows  a  mesh  laid  over  a  quarter  section  of  double-ridged  waveguide. 


Vertical  symmetry 


plane  (odd) 


Fig.  4.  Cross  section  of  a  double-ridged  waveguide. 
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Fig.  5.  Computer  mesh  for  a  quarter  section  of  a  double-ridged 

waveguide.  Darkened  lines  indicate  conducting  boundaries. 

The  nodes  of  the  mesh  are  described  to  the  computer  in  matrix  form. 
They  are  the  points  at  which  the  scalar  field  ^[10j  °  Hz  is  computed. 

The  field  form  of  the  Helmholtz  wave  equation  must  be  replaced  by 
a  discrete  form  in  accordance  with  the  physical  mesh.  This  is  done  by 
expanding  the  Laplacian  operator  in  Eq .  36  in  terms  of  a  truncated  two- 
dimensional  Taylor's  series.  The  result  for  node  5,  which  has  four 
nearest  neighbors  that  do  not  touch  the  conducting  boundaries  is 

♦j  +  +  ^3  +  \  +  f(uph)2  ~  ^5  =  0  (37) 

where  Up  is  an  approximate  numerical  eigenvalue  and  h  is  the  distance 
between  nodes.  The  form  of  Eq .  37  changes  for  nodes  on  conducting 

boundaries  and  symmetry  planes. 
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Symmetry  plane  and  boundary  points  are  handled  in  such  a  way  as  to 


j  enforce  two  rules:  1)  the  normal  derivative  boundary  condition  for  TE 

modes  and  2)  the  total  longitudinal  flux  of  the  guide  must  be  zero.  For 

t 

\  the  TEjq  mode,  the  normal  derivative  of  the  Hz  field  must  be  zero  along 

|  both  the  odd  symmetry  plane  and  the  waveguide  walls.  This  is  equivalent 

.  to  requiring  that 


o 

dn 


(38) 


since  Hz  is  proportional  to  iji.  Equation  38  can  be  written  in  its  cen¬ 
tral  finite  difference  form  for  node  U'g  as 


3*8  -  *7  -  *6 


3n 


2h 


=  0 


(39) 


which  means 

4^  =  4>7  (AO) 

The  node's  exterior  to  the  conducting  boundaries  and  odd  symmetry  plane 
make  it  possible  to  numerically  enforce  the  general  form  of  Eq .  AO. 
Quite  a  different  tactic  is  used  to  handle  the  even  symmetry  plane. 

In  order  to  prevent  the  numerical  method  from  giving  solutions 
representative  of  the  impossible  TEqq  mode,  the  total  longitudinal  flux 
of  the  guide  must  be  zero  (V  •  B  =  0).  This  condition  can  be  enforced 
by  requiring  that 
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at  the  even  symmetry  plane.  The  program  in  Appendix  B  encodes  the 
second  requirement  by  setting  the  average  value  of  the  interior  nodes  (h 
away  from  this  plane)  equal  to  0.5. 

By  writing  Eq .  37  at  the  m  (interior,  conducting  boundary,  even 
and  odd  symmetry)  nodes,  one  obtains  m  equations  involving  m  +  1 
unknowns.  This  is  now  a  matrix  eigenvalue  problem, 


Ai|i  * 


where  the  eigenvalues  are 


x  ■  (uph)2  =  (2lTh/Xc)2 


and  i  is  the  cutoff  wavelength.  The  eigenvector  will  be  made  up  of  the 
c 

field  points  <Ji  ,  4*~,  ...,  4*  •  The  matrix  eigenvalue  problem  can  now 

1  L  3  m 

be  solved  on  a  computer. 

The  matrix  eigenvalue  problem  is  solved  using  a  version  of  the 


inverse  power  method  called  doubly  iterative  successive  over-relaxation. 
In  this  technique,  computed  values  of  ^  are  used  to  obtain  an  approxima¬ 


tion  for  Up.  The  process  continues  until  an  iteration  is  reached  for 
which  the  previous  and  Up  are  within  some  user  specified  range  of  the 
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present  41  and  Up.  The  process  begins  with  a  guess  for  Up  that  is  pref¬ 
erably  less  than  the  actual  Up.  Equation  37  is  computed  at  each  node 
and  is  found  not  to  equal  zero.  Instead,  numbers  called  residuals  (R) 

are  obtained.  They  are  used  to  sequentially  replace  each  value  of 
with  the  old  value  plus  a  correction  dependent  on  the  residual, 


,  new 


.old 


wR 


(4  -  u2h2) 
^  P  ' 


(44) 


where  1  <  u>  <  2  is  the  over-relaxation  factor.  There  is  an  optimum 
value  of  U)  which  gives  a  final  solution  in  the  least  number  of  itera¬ 
tions.  Unfortunately,  it  must  be  obtained  empirically.  The  new  field 
values  are,  of  course,  wrong  since  a  wrong  initial  guess  of  Up  was  used. 

The  Rayleigh  coefficient  concept  uses  the  new  values  of  1 p  to 

2 

obtain  a  more  accurate  value  of  u_.  The  Rayleigh  coefficient  u  is 

P  p 

obtained  by  integrating  over  the  waveguide  cross  section  as 


2 

u 

P 


,  >2 

-/  iJA^da 
/  i|>2da 


(45) 


Note  that  the  discrete  form  of  da  is  different  for  the  nodes  A,  B,  C, 

and  D  shown  in  Fig.  5.  The  area  element  da  becomes  Aa  .  The  code  in 

2  2  2  2 

Appendix  B  assigns  Aa  values  of  0.25  h  ,  0.5  h  ,  0.75  h  ,  and  h  to 

nodes  like  A,  B,  C,  and  D,  respectively.  The  finite  difference  equiva¬ 


lent  of  Eq.  45  is 


The  summations  are  over  the  interior  and  boundary  points  of  the  guide 


cross  section.  The  to  's  are  implicitly  multiplied  by  the  appropriate 

i  *  J 

value  of  Aa .  Both  Fig.  6  and  the  procedure  outlined  below  describe  the 

2 

doubly  iterative  calculation  scheme  used  to  obtain  (u  h)  and  ip .  The 
final  eigenfunction  must  be  scaled  to  match  Solymar's  normalization  as 
shown  in  Appendix  E. 

2 

1.  Assume  initial  values  of  fu  h")  and 

p 

2 .  Use  Eqs.  37,  40,  41,  and  44  in  several  relaxation  passes  to 
relax  the  point  potential-function  values  to  a  reasonable 
degree . 

2 

3.  Use  Eq .  46  to  obtain  an  improved  estimate  of  (u^h)  . 

4.  All  nodes  exterior  to  the  conducting  boundaries  and  odd  sym¬ 

metry  plane  are  set  equal  to  the  interior  points  opposite 
them  (i.e.,  ^  Nodes  along  the  even  symmetry  plane  are 

held  at  ^  =  0,  while  nodes  1  mesh  unit  away  are  held  at  an 
average  value  of  0.5. 


5.  Iterations  will  be  stopped  when  both  the  largest  field  residue 


and  the  relative  difference  between  the  two  most  recent  values 
2 

of  (u  h)  are  less  than  their  convergence  criterion. 


Discretization  of  both  the  double-ridged  waveguide  cross  section 
and  the  Helmholtz  Wave  Equation  has  made  it  possible  to  describe  the 
problem  to  a  computer.  By  enforcing  the  appropriate  boundary  conditions 
and  applying  the  five-step  solution  procedure,  the  TEjq  mode  eigenvalue 
Ujq  and  eigenfunction  <Ji^g  can  be  obtained  at  any  particular  cross 
section  within  the  transition.  These  numbers  are  then  used  to  find  the 
coefficients  that  couple  the  incident  and  reflected  parts  of  the  TE^g 
mode . 


2 .  Computing  Coupling  Coefficients  for  the  Dominant  Mode 

The  problem  at  hand  is  to  find  the  coupling  coefficients  that  are 
needed  to  describe  TE^g  mode  propagation  in  a  double-ridged  waveguide. 
The  coupling  coefficients  8^  and  tc^  can  be  readily  computed  using  Eq . 
16  and  a  knowledge  of  u^g  (h^Qj  in  Eq.  16).  According  to  Eqs.  14  and 

A. 10,  S[io][iO]  is  the  only  coefficient  needed  since  [10]  iS 
zero.  Equation  A. 6  shows  that  this  coefficient  can  be  written  as  an 
integral  around  the  waveguide  boundary  C(x,y,z), 


S [ 10] [ 10] 


-  y  J  tan  0 

C(z) 


M  - 


(47) 


where  ds  is  an  element  of  length  along  C(x,y,z)  and  0  is  defined  in  Fig. 
1.  The  four  factors  which  contribute  to  a  successful  computation  of 
S[ 10] [ 10]  3re:  correctly  assigning  a  value  to  tan  0,  2)  finding 


tangential  derivatives  of  4^  at  the  boundaries,  3)  accounting  for  corners 
while  integrating  along  the  boundary  and  4)  choosing  an  appropriate 


value  for  the  node  spacing  h.  The  first  function  in  the  integrand  of 
Eq.  47  is  tan  0. 

a.  The  Definition  of  Tan  9 

The  value  of  tan  9  depends  upon  the  x-y  position  on  the  boundary 
and  the  axial  location  z  of  the  cross  section.  For  example,  the  program 
in  Appendix  B  requires  four  values  of  tan  0  to  describe  the  taper  flare, 
one  along  each  of  the  boundary  sections  shown  in  Fig.  7.  Equation  A. 11 
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Fig.  7.  Quarter  section  of  a  double-ridged  waveguide 
showing  the  _four  conducting  boundary  edges 
along  which  [10]  *S  coinPute^* 

defines  tan  0  using  the  notation  of  Fig.  1.  According  to  this  defini¬ 
tion,  tan  0  is  negative  at  points  where  the  boundary  slopes  toward  the 
z-axis.  Conversely,  tan  9  is  positive  for  those  points  where  the  bound¬ 
ary  slopes  away  from  the  axis.  The  second  function  in  the  integrand  of 


Eq.  47  depends  upon  <Ji. 


Cubic  splines  were  used  to  evaluate  the  tangential  derivative  of  4* 


along  the  four  boundary  lines  shown  in  Fig.  7.  In  this  method,  a  cubic 
polynomial  is  fit  to  the  boundary  field  data  set  [(Sj»4>j).  •••» 

(s  .*)]  as 

n  n 

*<s)  -  R(s)  -  *4  +  Bi  (s  -  S±)  +  CA  (s  -  Si)2  +  Di(s  -  S^3  (48) 

where  s  is  any  physical  point  along  C(x,y,z)  defined  on  the  interval 
between  and  S^.  The  spline  coefficients  B^,  and  are  computed 
from  the  data  set.  The  tangential  derivative  of  the  boundary 

field  can  be  obtained  for  any  boundary  point  s  by  evaluating  the  differ¬ 
entiated  form  of  Eq .  48, 

If :  ir1  •  Bi  *  2ct  *  h) +  3  Di  - s,)2  <4?> 

The  beginning  and  end  points  of  the  four  [s^,^]  data  sets  coincide  with 
the  end  points  of  the  boundary  sections  shown  in  Fig.  7.  This  segmenta¬ 
tion  of  9i|i/3s  was  necessary  since  tan  0  is  discontinuous  at  the  wave¬ 
guide  corners. 

c.  Dealing  with  Corners 

In  order  to  avoid  problems  with  a  discontinuous  integrand  at  the 
corners,  Eq .  47  was  split  into  four  parts.  The  Gauss  -luadrature  inte¬ 
gration  algorithm  was  used  along  each  of  the  four  boundary  segments. 


The  sum  of  these  integrals  was  then  multiplied  by  four  to  account  for 
the  entire  boundary.  Since  for  ridged  waveguides,  there  is  no  analyti¬ 
cal  solution  for  ®[io][10]  an<*  c^us  no  wa^  checking  computed  values, 
the  simpler  case  of  a  rectangular  waveguide  was  tested. 


d.  Dependency  of  _  j_0]  [10]  °n  ^ 

The  accuracy  of  the  ®[io] [10]  ca^cu^at^on  depends  upon  both  the 
precision  of  the  computer  used  and  the  node-to-node  spacing  h.  The  pro¬ 
gram  in  Appendix  B  was  used  to  compute  ^^q][10]  ^or  t*ie  rectanRular 
waveguide  shown  in  Fig.  8. 


y 


Fig.  8.  Rectangular  waveguide  used  to  test  the 
ridged  waveguide  program,  b/a  *  0.5. 


Table  1  shows  how  ^[iq][10]  aPProac^es  P*16 
(derived  in  Appendix  D)  of  0.1  as  the  ratio 


analytically  obtained 
of  h/a  decreases. 


value 
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Table  1.  Dependency  of  S 


on  h/a 


[10] [10] 


h/a 

s[ 10] [10] 

Error  (%) 

0.05 

0.0912 

8.8 

0.025 

0.0953 

4.7 

0.0125 

0.097589 

2.41 

0.00625 

0.098837 

1.16 

0.003125 

0.099640 

0.36 

0.0015625 

0.096163 

3.837 

The  optimum  accuracy  of  0.36  percent  error  is  a  result  of  machine 
precision.  The  results  in  Table  1  were  calculated  in  single  precision 
using  an  HP1000  computer.  Had  the  computations  been  performed  using 
double  precision,  the  optimum  value  of  h/a  would  have  been  much  smaller 
and  the  accuracy  of  ^[iq] [10]  muc^  8reater* 

The  results  of  this  computer  program  test  can  be  extended  to  the 
ridged  waveguide  case.  The  finite  difference  numerical  algorithm  used 
to  obtain  the  eigenvalues  and  eigenvectors  has  been  successfuly  tested 
against  Cohn's9  results  for  ridged  waveguides.  Since  the  computer 
program's  accuracy  is  limited  by  the  working  precision  of  the  machine, 
it  is  reasonable  to  assume  that  results  accurate  to  within  0.3  percent 
can  be  obtained  for  waveguides  containing  ridges. 

The  TEjO  mode  coupling  coefficients  k,  6,  and  S[^o][10]  *iave  been 
shown  to  depend  upon  the  geometrical  and  electrical  characteristics  of 
the  particular  cross  section  in  question.  These  coefficients  are  now 
one  step  away  from  their  final  form. 
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Piecewise  Continuous  Coupling  Coefficients 


In  their  final  form,  the  coupling  coefficients  are  represented  by 
piecewise  continuous  functions  of  axial  position  z.  This  is  accom¬ 
plished  by  computing  them  at  discrete  points  Z^  along  the  transition. 

The  data  sets  [Ki*zi]  an^  [ S [  1 0]  [  10] 1  ,Z1  ^  are  t0  a  cubic 

spline  similar  to  Eq.  48.  The  number  of  points  chosen  to  represent  a 
specific  transition  is  left  to  one's  discretion.  Large  changes  in 
waveguide  geometry  which  occur  within  one  axial  wavelength  will  necessi¬ 
tate  a  finer  discretization  in  order  to  accurately  capture  the 
behavior  of  the  coefficients. 

In  summary,  the  numerical  design  tool  for  the  TE^q  mode  double- 
ridged  waveguide  has  been  developed  to  the  point  of  representing  the 
coupling  coefficients  as  piecewise  continuous  cubic  splines.  The  normal 
mode  equations  can  now  be  solved  for  the  transition  scattering  matrix  S. 

B.  TEjq  Scattering  Matrix  for  a  Double-Ridged  Taper 

According  to  Section  1,  the  TE^q  mode  system  of  coupled  differen¬ 
tial  equations  can  be  transformed  into  the  transition  scattering  matrix 
in  three  steps:  1)  convert  complex  equations  to  real  ones,  2)  solve  the 
real  equations  to  obtain  the  transmission  matrix  T  and  3)  use  T  to 

n 

obtain  the  scattering  matrix  S. 

1.  TE^q  Mode  Equation  Conversion  Complex  to  Real 

In  order  to  use  the  differential  equation  solver  DESOLV10  listed 
in  Appendix  B,  the  TEjq  form  of  Eq .  14, 


1  d(lntc) 

2~dT~A  +  S 


iiOi£iA%s-  A+ 


must  be  converted  into  an  equivalent  real  matrix  equation.  The  TEjq 
mode  bracket  notation  [10]  used  in  Appendix  A  has  been  dropped  for  the 
sake  of  clarity.  Equation  50  may  be  expressed  in  the  matrix  notation  of 
Eq.  17  as 


i.  ln<)‘ 
2  dz 


JL  d(ln<) 

2  d2 


where  the  +  and  -  superscripts  denote  forward  and  backward  waves, 


respectively.  To  simplify  the  notation,  Eq .  51  can  be  converted  into 
the  following  form 


Since  both  the  A's  and  the  M's  have  real  and  imaginary  parts,  they 
may  be  rewritten  as 


-  3b  - 


-'W*'  ■■'‘/-'S'. % ■"/ :*v 


-.w  •* 


A±r+  j 


,±1 


M  =  Mr  +  jM1  (53) 

mn  mn  mn 


where  the  r  and  i  superscripts  refer  to  real  and  imaginary  parts,  and 
the  m  and  n  subscripts  refer  to  elements  in  the  m1"^1  row  and  nc^  col¬ 
umn.  Applying  Eq .  53  to  row  2  of  Eq.  52  gives 


-r 


+  j  A 


-i 


)  =  (Mjj  +  J  H 


21 


)(A+r  +  j 


A+1) 


(^22  +  J  r  +  j  A 


-i 


‘22 


(54) 


Carrying  out  the  multiplication  and  equating  real  and  imaginary  parts 
yields 


A  r  =  A  -  M2 j  A  +  M22  A  -  M22  A 

A-1  =  A+r  +  HT2l  A+i  +  A_r  +  M^2  A_i  (55) 


Performing  this  set  of  operations  on  row  1  of  Eq .  52  will  give  a  set  of 
equations  similar  to  Eq .  55  with  M2  ^  and  M22  replaced  by  Mjj  and  M^> 
respectively: 


d_ 

dz 


+r 


d_ 

dz 


+i 


A+r  -  M1 
A  Hjj 


^  *  Hj, 


a+1  +  m[2  A_r  -  n}2  a'1 
A+‘  ♦  M;2  A'r  ♦  «;2  A-‘ 


(56) 
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Substituting  Eq .  58  into  Eq .  57  results  in  the  real  matrix  form  of  the 


TEifl  mode  coupled  differential  equations. 


,  1  d ( Inic(z ) 

S  <2)-  2  - 


c;-/2n_  I  dlln£C2J 

b  KZ)  2  dz 


-6(z) 


S  (z)  - 


_1  d  ( ln<(  z  )  ) 
2  dz 


+  6(z) 


Equation  59  can  now  be  used  to  obtain  the  transmission  matrix. 


2 .  Formulating  the  Transmission  Matrix 


Extracting  the  TE,n  mode  transmission  matrix  from  Eq .  59  is  a  two- 


step  process:  1)  solve  it  twice  using  orthogonal  mode  amplitude  initial 


condition  vectors  and  2)  express  the  transition's  mode  amplitudes  at  the 


output  in  terms  of  those  at  the  input. 


The  orthogonal  initial  conditions  y^  and  shown  in  Eq .  60  repre¬ 


sent  forward  and  backward  waves  at  the  transition  input  with  unit  ampli¬ 


tude  and  zero  phase. 


^  ^  „  av  ;X;£ 


yj(o)  = 


y2(0)  = 


The  first  two  rows  represent  the  real  and  imaginary  parts  of  the  forward 


wave;  likewise,  the  second  two  rows  represent  the  backward  wave.  Solv¬ 


ing  Eq.  59  with  these  initial  conditions  yields  the  following  linearly 


independent  solution  vectors  at  z  =  L, 


yjCU 


y2(L) 


These  initial  condition  and  solution  vectors  are  algebraically 


transformed  into  the  transmission  matrix  as  follows.  The  real  form  of 


the  problem  is  returned  to  its  original  complex  form  by  rewriting  Eqs. 


60  and  61  in  terms  of  the  complex  variables  u,  T  and  x. 


iXjlA  aJI  kA  -Ml 


(L)  =  l  C  x  (L)  =  Cx 


11 

T 

12 

21 

T 

22 

=  T  C 


The  fact  that  C  is  common  to  both  Eqs.  63  and  64,  makes  it  possible  to 
express  forward  and  backward  mode  amplitudes  at  z  =  L  in  terms  of  those 
at  z  =  0. 


x  (L)  =  T  x  (0) 


Equation  65  is  in  the  form  of  Eq.  31  where  T  is  the  TE^q  mode  transmis¬ 
sion  matrix.  The  transmission  matrix  is  finally  used  to  obtain  S. 


3.  Transmission  Matrix  to  Scattering  Matrix 


The  TEjq  mode  scattering  matrix  S  is  obtained  by  rearranging  Eq . 
65  in  terms  of  the  scattering  notation  of  incident  and  reflected 
waves.  The  process  begins  by  rewriting  x  in  Eq .  65  in  terms  of  A+  and 

A", 


A+(L) 


t12  a+(0) 


A  (L) 


T22  J  LA  (0) 


n  .* 

*■.  .  • 
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Figure  9  illustrates  Eq .  66  in  terms  of  the  scattering  notation  of  Fig. 


A  (0)  =  al 


A  (L)  =  b2 


Tapered  transition 


A” (0)  =  bl 


A" (L)  =  a2 


Fig.  9.  A  single  mode  illustration  expressing  A+  -  A  nota¬ 
tion  in  terms  of  incident  and  reflected  waves. 


Equation  66  can  be  rewritten  in  scattering  notation  as 


Notice  that  Eq .  67  is  in  the  form  of  Eq .  32,  and,  as  expecteu  for  a 


single  mode  analysis,  the  matrix  and  vector  notations  are  gone.  The 


general  expression  for  transforming  transition  matrices  into  their 


corresponding  scattering  matrix  given  by  Eq .  34  can  be  applied  to  Eq .  67 


with  the  following  result, 


.  - 
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T  -  T  T 
11  12  22  21 


T  T  * 
12  22 


a  1 


a2 


(68) 


This  may  be  rewritten  in  the  form  of  Eq .  35  as 


>  =  > 
b  =  S  a 
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spline. 

S  is  the  TE^q  mode  scattering  matrix  of  an  arbitrarily  tapered  double- 
ridged  waveguide  transition  having  cross  sections  with  quarter-waveguide 
symmetry.  The  program  RIVSWR  in  Appendix  B  has  been  designed  to  imple¬ 
ment  this  single-mode  version  of  the  multimode  analysis  technique. 

To  summarize,  two  major  aspects  of  numerically  obtaining  a  scat¬ 
tering  matrix  have  been  presented.  First,  in  order  to  solve  the  TE^q 
mode  coupled  system  of  differential  equations,  the  coupling  coeffici- 

had  to  be  known  as  continuous  functions  of 


f 10] [10J 


using  the  eigenvalue  and  eigenfunction  of  the  TE10  mode.  The  finite- 
difference  inverse  iterative  power  method  was  used  to  solve  the  matrix- 
eigenvalue  problem.  Gaussian  integration  was  used  to  obtain  S. 


[10] [10] 

from  the  waveguide  boundary  fields.  Second,  with  the  coupling  coeffici¬ 


ents  in  hand,  the  routine  DESOLV  was  used  with  mutually  orthogonal 
initial  condition  vectors  to  find  linearly  independent  solution  vectors 
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IV.  EXPERIMENTAL  VERIFICATION  OF  PROGRAM  VSWR  PREDICTIONS 


Comparisons  between  measured  and  computed  results  show  that  the 
taper  analysis  technique  presented  herein  can  be  used  to  accurately 
predict  transition  performance.  In  order  to  use  terms  which  better  suit 
measured  data,  this  section  places  emphasis  on  VSWR  (computed  from  S- 
parameters).  Comparisons  are  made  between  computed  and  measured  VSWR 
versus  frequency  profiles  for  two  linearly  tapered  unridged  transi¬ 
tions.  The  comparisons  show  that  the  code  is  valid  for  these  geome¬ 
tries.  A  detailed  explanation  is  given  regarding  how  the  code  was  used 
to  model  a  cosine  impedance  transition  tapering  from  rectangular  to 
double-ridged  waveguide.  Measurements  made  on  a  cosine  impedance  taper 
show  that  the  code  accurately  models  double-ridged  transitions  with 
nonlinear  tapers. 


A.  Two  Linearly  Tapered  Transitions  in  Rectangular  Waveguide 

The  work  of  S.  S.  Saad11  and  Z.  Wenxin12  is  compared  to  results 
generated  by  the  program  in  Appendix  B;  within  experimental  error,  the 
code  accurately  models  unridged  transition  performance.  The  VSWR  pro¬ 
file  reported  by  Saad  for  height  tapered  transitions  agreed  with  the 
code's  predictions.  Similarly,  the  code  accurately  predicted  Wenxin's 


VSWR  profile  for  a  transition  linearly  tapered  in  both  height  and  width. 


1 .  Computed  Versus  Measured:  Normal  Mode,  Saad  and  Young 


According  to  the  VSWR  data  computed  by  Saad  and  measured  by  L. 
Young,13  the  code  accurately  models  dominant  mode  behavior  in  linearly 
height  tapered  rectangular  waveguides.  Figure  10  shows  the  symmetrical 


Fig.  10.  Symmetrical  linear  height  tapered 
transition  in  WR-650  waveguide. 

linear  taper  originally  analyzed  by  Young.  He  measured  the  taper's  VSWR 
by  placing  it  back  to  back  with  a  quarter  wave  transformer  having  a  1.01 
maximum  VSWR.  His  measurements  are  shown  in  Fig.  11,  along  with  Saad’s 
numerical  solution  and  the  code's  predictions. 

As  these  VSWR  profiles  show,  the  code  is  capable  of  modeling 
tapers  like  the  one  shown  in  Fig.  10.  Considering  possible  differences 
between  computer  programs  (precision,  algorithms,  error  tolerances, 
etc.),  the  VSWR  profiles  computed  by  Saad  and  the  normal  mode  code  are 
in  good  agreement.  The  code  is  not  a  complete  model  of  the  taper; 
losses,  higher  order  modes,  and  mechanical  imperfections  are  not  taken 
into  account.  Likewise,  the  measured  data  are  not  error  free.  With 
these  facts  in  mind,  the  agreement  between  measured  and  computed  VSWR 
profiles  is  quite  acceptable.  Results  similar  to  these  have  also  been 
obtained  for  doubly  tapered  rectangular  waveguides. 


Computed  and  measured  VSWR  profiles  for  the  transition  shown  in  Fit*.  10 


2.  Computed  Versus  Measured:  Normal  Mode,  Wenxln  and  Johnson 


The  VSWR  data  measured  by  Johnson*4  and  computed  by  Wenxln  show 
that  the  normal  mode  code  is  also  capable  of  modeling  the  dominant  mode 
performance  of  unridged  transitions  whose  broad  and  narrow  sides  are 
linearly  tapered.  In  1959,  Johnson  measured  the  VSWR  of  the  doubly- 
tapered  transition  shown  in  Fig.  12. 


Fig.  12.  Doubly-tapered  rectangular  waveguide 
analyzed  by  Johnson  and  Wenxin. 

The  measured  data  for  this  transition  show  that  the  code's  VSWR 
predictions  will  be  low  for  frequencies  above  which  higher  order  mode 
propagation  occurs.  The  measured  and  computed  VSWR  profiles  are  shown 
in  Fig.  13.  Notice  the  difference  between  the  predictions  of  Wenxin, 
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the  normal  mode  code  and  measured  data  for  frequencies  above  about  9.8 


GHz.  The  predicted  VSWR  is  low  for  this  portion  of  the  curve.  The  TEqj 
mode  becomes  transmissible  within  the  taper  at  9.8  GHz.  Since  its 
effect  on  the  TEjq  mode  is  not  included  in  the  numerical  model,  the 
theoretical  prediction  of  VSWR  should  be  lower  than  the  measured  one. 
With  the  exception  of  Schelkunoff ,  1  ^  the  computed  VSWR  profiles  were 
very  accurate  below  9.8  GHz.  Results  similar  to  those  presented  for 
unridged  waveguide  transitions  have  also  been  obtained  for  ridged  ones. 

B.  A  Cosine  Impedance  Transition  in  Double-Ridged  Waveguide 

This  work  culminates  in  the  ensuing  paragraphs  where  the  agreement 
between  theory  and  experiment  shows  that  the  normal  mode  technique  is 
capable  of  successfully  predicting  the  VSWR  profiles"  of  nonlinear  wave¬ 
guide  tapers.  A  detailed  example  is  given  of  how  the  normal  mode  code 
RIVSWR  (Appendix  B)  was  used  to  transform  the  physical  dimensions  of  a 
cosine  impedance  transition  (WR-90  to  WRD-750)  into  a  VSWR  versus  fre¬ 
quency  profile  for  the  dominant  mode.  Network  analysis,  time  domain 
ref lectometry  and  inverse  Fourier  transforms  are  used  to  obtain  measured 
data  that  compare  well  with  the  code’s  prediction. 

1.  Transforming  Waveguide  Dimensions  into  a  VSWR  Profile 

In  order  to  run  the  code,  the  user  must  create  a  data  file  which 
provides  an  accurate  discretized  description  of  the  taper's  boundary 
(RSIZ.DAT) .  The  following  example  shows  1)  how  this  file  was  created 
for  the  cosine  taper  and  2)  a  sample  run  with  the  resulting  VSWR  pro¬ 
file. 


A  cosine  impedance  function  was  chosen  for  this  example  since  it 


can  be  used  to  make  short  low  VSWR  transitions.  The  function  is  given 


Zq(z)  =  (Z1Z2)1/2  exp  -  y  In  [{z^/Z^)\  cos  (nz/L) 


Z^  and  Z2  are  the  respective  characteristic  impedances  of  the  WRD-750 


and  WR-90  ends  of  the  taper  with  length  L  =  1  inch.  A  plot  of  Eq .  69  is 


shown  in  Fig.  14.  A  definition  of  impedance  in  terms  of  waveguide 


dimensions  was  used  to  impose  this  profile  upon  the  transition, 


Hoefer's16  voltage  to  current  based  definition  of  ridged  waveguide 


impedance  was  used  to  find  an  axial  profile  for  ridge  height.  Figure  15 


shows  the  notation  Hoefer  used  to  define  the  impedance 


Z  -  Z  1  -  f X/X  V 
o  o00  v  cr' 


where 


120*  (b/*cr) 


b  v  *s 

d  510  ~  + 
cr 


o  .  .  *  b  a  -  s 

—  +  tan  y  —  — T" 

o  cr  \ 


b  ,  ,  4  /  n  -  f  b  V  '  b  - 
TCa  '-'tt  1+7  1  +  °-2  r^rlncsc2b- 
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Fig.  15.  Dimensions  of  a  double-ridged 

waveguide  as  defined  by  Hoefer. 

b/X  is  the  normalized  cutoff  frequency.  The  normalized  susceptance  is 
cr 

approximately 


Bq/Yo  -  (2b/Xcr)  In  esc  ~  (73) 

In  order  to  simultaneously  solve  Eqs.  69  and  70  for  d  at  a  number  of 
axial  positions,  s  was  kept  constant  (0.73”)  and  a  and  b  were  linearly 
tapered  from  one  end  to  the  other.  With  a,  b,  s  and  ZQ  specified  at 
0.01  inch  intervals  in  axial  position,  a  root  finding  routine  was  used 
to  solve  Eq .  70  for  d  at  101  points  along  the  taper.  A  profile  of  the 
ridge  height  was  obtained  from  d  and  is  shown  in  Fig.  16.  In  addition 
to  specifying  the  geometry  of  the  transition,  RSIZ.DAT  must  contain 
information  about  the  slopes  of  the  waveguide  boundaries. 

Unlike  the  other  waveguide  boundaries,  a  least  squares  fit  was 


applied  to  the  ridge  height  profile  in  order  to  obtain  a  smooth  slope. 


The  slope  of  the  taper  in  b  from  the  input  bl  to  the  output  b2  was 
obtained  as 


tan  02  = 


b2  -  bl 


2L 


0.321  -  0.4 


2(1) 


=  -  0.0395 


(74) 


Similarly,  the  slope  of  the  taper  in  a  (tan  01)  was  calculated  to  be 
-0.1045.  Since  the  ridge  width  was  held  constant,  the  slope  in 
s  (tan  03)  is  zero. 

An  eighth  order  fit  on  the  computed  boundary  data  for  ridge  height 
was  used  to  obtain  its  slope  as  a  function  of  axial  position.  The 
numerical  inaccuracies  of  the  root  finding  computations  were  smoothed 
away  by  the  least  squares  fit.  The  fit  gives  the  ridge  position  h'  with 
respect  to  the  waveguide  axis  as  shown  in  Fig.  17.  The  fit  function  is 
given  by 

h ' ( z )  «  0. 138658z  +  0.408664 z2  -  1.02566z3  +  0.897162z4 

-  0.225208z5  +  0.022618z6  -  0.228794z7  +  0.144703z8  (75) 


The  derivative  of  this  curve  describes  how  the  ridge  moves  away  from  the 
z'  line.  Its  negataive  is  the  slope  of  the  ridge  boundary  d  with 
respect  to  the  z-axis. 
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Fig.  17.  Ridge  height  with  respect  tc  a  line 
parallel  with  the  waveguide  axis. 


tan  04  =  -(0.138658  +  0.817328x  -  3.07698x2  +  3.58865xJ 


-  1.12604x4  +  0.135708x5  -  1.60156x6  +  1.15762x7)  (76) 


In  summary,  nine  data  points  are  needed  to  describe  the  waveguide 
boundary  at  an  axial  position;  a,  b,  d,  s,  z,  tan  01,  tan  02,  tan  03  and 
tan  04.  The  data  file  developed  for  the  WR-90  to  WRD-/50  transition  is 
shown  in  Appendix  C.  The  first  line  contains  the  number  of  axial  posi¬ 
tions  for  which  data  are  given.  Every  two  lines  thereafter  contain  the 
dimensions  and  tangent  data,  respectively.  This  file  was  used  by  RIVSWR 
(Appendix  B)  to  obtain  the  VSWR  profile. 


The  following  sample  run  of  RIVSWR  shows  how  to  input  data  and 
where  to  find  computed  results.  The  code  assumes  that  RSIZ.DAT  contains 
the  appropriate  data.  The  user  types  in  the  underlined  portions. 

$  RUN  RIVSWR 

ENTER  REFLECTION  COEFFICIENT  OF  SOURCE 

C0.,0.) 

ENTER  REFLECTION  COEFFICIENT  OF  LOAD 

(Q.,0.) 

ENTER  LOW  AND  HIGH  EDGES  OF  SWEEP  BAND  (GHz) 

8.4,18.0 

ENTER  #  OF  FREQUENCY  STEPS 
100 

EIGENVALUES  IN  EIGDAT.DAT?  TYPE  "1”  IF  SO 

2_ 

HOW  GOOD  SHOULD  THE  FIT  BE?  (INCHES) 

.001 

ERROR  OF  FIT  =  6.0239E-4  (INCHES)  H  =  2.1708E-3  (INCHES) 
ACCELERATION  FACTOR  W  <=  1  .9385 
CUTOFF  FREQUENCY  =  6.536667  GHz 

In  addition  to  the  users  guide  in  Appendix  B,  a  brief  explanation 
will  be  made  regarding  the  above  run.  If  the  user  wishes  to  model 
transition  performance  in  the  presence  of  load  and  source  mismatches, 
complex  reflection  coefficients  other  than  those  shown  may  be  entered. 
For  the  above  example,  the  code  will  attempt  to  fit  its  mesh  (which 
represents  a  quarter  of  the  waveguide)  to  within  0.001  inches  of  the 
waveguide  boundary.  This  represents  a  maximum  total  fit  error  of  0.002 
inches.  The  error  of  fit  is  limited  only  by  the  size  of  the  matrix  HZ 
in  RIVSWR.  The  printout  sequence  from  ERROR  OF  FIT  to  CUTOFF  FREQUENCY 
continues  until  all  the  cross  sections  of  RSIZ.DAT  have  been  analyzed. 
The  code  then  writes  the  frequency,  S-parameter  and  VSWR  data  to  the 
files  SPARAF1.DAT  and  PVSWR.DAT.  Figure  18  shows  the  VSWR  versus  fre¬ 
quency  profile  for  the  WR-90  to  WRD-750  taper.  As  will  be  seen  in  the 
following  pages,  this  computed  profile  agrees  well  with  the  measured 


Cosine  Taper  VSWR  Measurements 


The  capabilities  of  the  Hewlett  Packard  HP8510A  network  analyzer 
were  used  to  obtain  the  VSWR  profile  of  the  cosine  impedance  -  taper .  In 
addition  to  the  HP85lOA's  waveguide  calibration  kit,  its  time  domain 
ref lectometry  and  inverse  Fourier  transform  functions  helped  make  accu¬ 
rate  VSWR  measurements  of  the  cosine  impedance  taper.  Figure  19  shows 
two  views  of  the  electrof ormed  taper. 


a.  The  Experimental  Setut 


A  WRD-750  sliding  load  and  offset  shorts  were  used  to  calibrate 
the  system  out  to  the  test  plane.  The  test  plane  was  the  open  end  of  a 
WRD-750  waveguide.  The  other  end  of  the  waveguide  was  attached  to  the 
system  by  a  coax  to  waveguide  transition.  The  phase  and  nu^rnitude  of 
Sjj  was  measured  by  the  system  for  two  waveguide  shorts  (0.256  and  0.768 
inch  offsets)  from  8.4  to  12.4  GHz.  Sjj  for  the  sliding  load  was  mea¬ 
sured  across  the  same  band  for  several  different  load  locations.  This 
procedure  was  repeated  over  the  12.4  to  18  GHz  band  with  0.118  and  0.354 
inch  offset  shorts.  Figure  20  shows  a  block  diagram  of  the  experimental 
set-up.  The  load  was  separated  from  the  cosine  taper  by  18  inches  so 
that  reflections  from  each  would  be  well  separated  in  time. 
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Time  domain  ref lectometry  and  inverse  Fourier  transforms  were 
successfully  used  to  filter  out  the  load's  effect  on  VSWR.  '  The  time 
domain  ref  lectometry  data  taken  over  the  low  (8.4-12.4  GHz)  band  is 
shown  in  Fig.  21.  The  highest  peak  (first  in  time)  corresponds  to  the 
left  edge  of  the  taper  (WRD-750)  and  the  second  peak  corresponds  to  the 
right  edge  (WR-90).  As  expected,  the  load  response  (third  peak)  is  very 
distinct  from  the  other  two.  A  similar  response  was  obtained  for  the 
high  band  (12.4-18  GHz).  The  VSWR  of  the  taper  with  an  ideal  load  can 
now  be  approximated  by  neglecting  the  load  response. 

By  taking  the  inverse  Fourier  transform  of  a  gated  portion  of  the 
time  response  curve,  the  load's  effect  on  measured  taper  VSWR  was  elimi¬ 
nated.  This  is  clearly  shown  in  Fig.  22  where  the  rippled  and  smooth 
curves  correspond  to  measured  and  modified  VSWR  data,  respectively.  The 
smooth  curve  (8.4-12.4  GHz)  corresponds  to  the  gated  portion  (between 
markers)  of  the  time  response  curve  in  Fig.  21.  The  smooth  portion  of 
the  high  band  curve  was  obtained  in  the  same  manner.  The  smooth  curve 
is  used  to  represent  measured  data  in  the  comparison  with  computed 
predictions . 


3.  Computed  Versus  Measured  VSWR 


The  good  agreement  between  measured  and  computed  VSWR  provides  the 
final  piece  of  evidence  in  support  of  RIVSWR  and  the  normal  mode  tech¬ 


nique.  Figure  23  shows  these  profiles  for  the  cosine  impedance  taper. 


es 


Three  reasons  for  the  difference  between  these  profiles  are  1) 
physical  modeling  error,  2)  higher  order  mode  propagation  and  3)  mea¬ 
surement  error.  In  the  sample  run  of  RIVSWR,  the  code  was  forced  to  fit 
each  waveguide  cross  section  to  within  -0.002  inches.  On  the  average, 
the  dimensions  of  the  taper  deviated  from  design  by  -0.001  inches. 
These  physical  modeling  errors  accompanied  by  the  numerical  errors 
previously  mentioned  are  one  source  of  the  discrepancy  between  theory 
and  experiment. 

Higher  order  modes  are  partially  responsible  for  high  values  of 
measured  VSWR.  The  TE^  an<^  ^01  moc^es  become  transmissible  within  the 
transition  below  13.1  and  14.7  GHz,  respectively.  Since  RIVSWR  models 
dominate  mode  (TE  )  behavior,  its  VSWR  predictions  do  not  account  for 
the  effects  of  these  modes.  Below  about  11  GHz,  these  effects  are  not 
present  and  measurement  errors  are  more  readily  identified. 

Since  the  HP8510A's  resolution  is  very  good,  it  is  difficult  to 
ascribe  differences  between  theory  and  experiment  to  the  measurement 
process.  The  HP8510A  used  calibration  data  and  internal  error  correct¬ 
ing  routines  to  obtain  a  resolution  of  about  43  dB.  This  represents  a 
possible  error  in  VSWR  of  about  0.013.  In  light  of  this  fact,  the 
maximum  deviation  between  theory  and  experiment  (from  8.4  to  12.4  GHz) 
is  about  0.05  in  VSWR.  From  a  practical  viewpoint,  this  error  and  the 
0.09  VSWR  error  at  18  GHz  is  quite  acceptable. 

To  summarize,  experimental  and  computed  VSWR  profiles  have  been 
presented  for  three  waveguide  tapers.  They  represent  the  ridged  and 
unridged,  linear  and  nonlinear  types  of  tapers.  In  all  three  cases, 
RIVSWR  accurately  predicted  the  measured  VSWR.  Aside  from  numerical. 


tolerance  and  measurement  error,  the  largest  errors  were  observed  for 
tapers  which  propagated  higher  order  modes.  The  twofold  purpose  of  this 
section  has  been  fulfilled.  First,  the  normal  mode  technique  presented 
in  previous  sections  has  been  shown  to  work.  Second,  the  example  of 
creating  RSIZ.DAT  and  running  the  code  for  a  nontrivially  tapered 
double-ridged  waveguide  makes  it  possible  for  other  workers  to  use 
RIVSWR  for  similar  applications. 


V.  CONCLUSION 


The  theoretical  and  numerical  aspects  of  developing  a  waveguide 
transition  design  tool  have  been  presented.  Normal  mode  and  finite 
difference  concepts  provided  a  foundation  for  the  code  RIVSWR.  It 
satisfactorily  predicts  the  dominant  mode  VSWR  of  tapers  in  rectangular 
and  double-ridged  waveguide. 

Calculated  and  measured  VSWRs  of  three  tapers  were  compared.  The 
first  two  were  linear  tapers  in  rectangular  waveguide.  The  third  was  a 
cosine  impedance  taper  in  double-ridged  waveguide.  In  each  case,  the 
code  predicted  VSWR  profiles  that  were  typically  within  5  percent  of 
measured  ones.  Within  the  operating  band  of  the  dominant  mode,  the 
error  was  mainly  attributed  to  tolerance  and  measurement  errors.  Above 
this  band,  higher  order  modes  propagate.  Since  the  code  does  not 
account  for  the  affect  these  modes  have  on  the  dominant  mode,  its  VSWR 
predictions  were  low  in  this  region. 

As  the  previous  paragraph  suggests,  the  code  could  be  improved  by 
modeling  higher  order  modes.  Additionally,  modes  passing  through  cutoff 
within  the  transition  could  be  included.  These  improvements  would  make 
it  possible  to  accurately  predict  the  VSWR  of  double-ridged  waveguides 


having  wide  flare  angles. 


APPENDIX  A 


COUPLING  COEFFICIENTS 


b'. 


The  coupling  coefficients  for  the  V-I  and  A+  -  A-  formulations  are 
listed  as  T  and  S,5  respectively. 


T(i)(p) 


h(p)  9*"' 


2  2 

h, . .  -  h,  . 

(i)  (p) 


^  tan  9  - ds  ;  h  .  *  h 

C(z)  Sn  3n  (1) 


<p) 


(A.l) 


T(i)[p]  "  ° 


(A. 2) 


[i ] [p] 


ILL 


2  2  ^  tan  8  *m 

hf  ,  -  hf C( z )  llJ 


3% 


Ip]  ti] 


3n 


L2idS;  Nil  *  h|p] 


(A. 3) 


34'r,1 

r  r  .  ,  ,  .  =  -  <j>  tan  9  p 

111<p)  C(z) 


Ti - ST 


ds 


(A. 4) 


*  K«z> tan  9  ffi2 


( i )  (  i ) 


S(i)(i) 


ds 


(A. 5) 


C(z) 


^  r  ,•  l  r  z  l  =  S  r  .  ,  r  .  ,  =  “  t  ^  tan  ® 


ds 


(A. 6) 


APPENDIX  B 

DOMINANT  MODE  DOUBLE-RIDGED  WAVEGUIDE  PROGRAM  DESCRIPTION, 

FLOW  CHART  AND  FORTRAN  LISTING 

I.  PURPOSE 

This  program  computes  the  scattering  matrix  and  VSWR  of  a  double 
ridged  waveguide  taper  at  equally  spaced  intervals  in  the  designated 
frequency  band. 

II.  SCOPE 

The  program  implements  a  dominant  mode  version  of  the  A+  -  A- 
formulation  and  will  predict  low  VSWR  profiles  for  tapers  that  excite 
"strong"  higher  order  propagating  or  evanescent  modes. 

III.  METHOD 

The  finite  difference  method  is  used  to  compute  cross-section 
eigenvalues  and  eigenvectors.  The  Gauss  integration  formula  is  used  to 
evaluate  Solynar's  coupling  coefficient  J>[iq][io]*  Shampine’s10  sub¬ 
routine  DESOLV  integrates  the  normal  mode  eauations  using  the  modified 
divided  difference  form  of  the  Adams  Pece  formulas.  Finally,  the  steps 
outlined  in  Sections  II  and  III  are  implemented  to  obtain  the  scattering 
matrix  and  VSWR . 

IV.  ORGANIZATION 

This  Appendix  contains  a  description,  flow  chart  and  Fortran 
listings  of  the  RIVSWR  and  DESOLV.  R1VSWR  is  composed  of  three  major 
parts.  The  first  part  assumes  the  user  has  (1)  divided  the  transition 


into  a  representative  set  of  cross  sections  and  (2)  sequentially 
provided  the  standard  dimensions  of  each  cross  section  in  a  data  file 
(RSIZ.DAT).  This  data  is  used  to  compute  coupling  coefficient 
information  at  each  cross  section  and  is  stored  in  a  data  file 
(EIGDAT.DAT).  The  second  part  of  RIVSWR  uses  the  coefficient  data  and 


DESOLV  to  compute  the  scattering  matrix  and  VSWR  at  the  frequencies 
specified.  The  input  VSWR  (VSWR1 )-f requency  set  is  written  to  PVSWR.DAT 


at  each  frequency.  The  S-parameter  and  frequency  data  are  written  to 
SPARAM.DAT.  The  third  part  of  the  program  contains  subroutines  for 
integration,  differentiation,  cubic  spline  fitting  and  mesh  fitting. 


V.  USAGE 


The  important  aspects  of  using  this  code  fall  under  the  following 
categories:  input  variables,  output  variables,  accuracy  specifiers  and 
error  detectors. 


Input  Variables 

These  are  read  in  from  both  the  te.  ainal  and  the  RSIZ.DAT  or 
EIGDAT.DAT  data  files.  The  variables  read  from  the  terminal 
are  defined  below  in  the  order  one  would  enter  them. 
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Read  In  from  terminal 

a.  GAMAS  —  a  complex  variable  (real,  imaginary)  that 

expresses  the  reflection  coefficient  of  the 
source.  Enter  (0.,0.)  if  source  has  no  reflec¬ 
tion  coefficient. 

b.  GAMAL  —  a  complex  variable  that  expresses  the  reflection 

coefficient  of  the  load. 

c.  FMIN  & 

FtlAX  —  lower  and  upper  edges  of  the  frequency  band  over 
which  S-parameter  and  VSWR  information  are 
desired.  For  example,  if  the  band  is  2-10  GHz, 
then  FMIN  =  2  and  FMAX  -  10. 

d.  NFS  —  number  of  frequency  steps  the  band  is  to  be 

divided  into.  Note,  if  NFS  -  8  for  the  2-10  GHz 
band,  there  will  be  9  frequency  data  points. 

e.  I YES  —  is  1  if  EIGDAT.DAT  contains  information  about  the 

current  geometry,  otherwise  IYES  *  2. 

f.  RFIT  —  required  fit  (inches).  Applies  to  1/4  of  the 

waveguide;  the  fit  error  for  the  entire  waveguide 
will  be  less  than  2*RFIT. 

Read  in  from  RSIZ.DAT 

Appendix  C  shows  a  data  file  for  the  cosine  taper.  The 
first  row  specifies  the  number  of  cross  sections.  Every 
two  rows  thereafter  contain  information  about  the  dimen¬ 
sions  and  tan  0  of  a  particular  cross  section. 


AA  —  interior  width  of  standard  ridged  waveguide 

(inches) 

BB  —  interior  height  of  standard  ridged  waveguide 

(inches ) 

DD  —  distance  between  ridges  (inches) 

SS  —  ridge  width  (inches) 

ZLOC(KKK)  —  distance  the  KKKth  cross  section  is  from  the 
beginning  of  the  transition  (inches) 

Row  3 

TANTi(KKK)  —  is  the  tangent  of  the  taper  angle  on  the  ith 
side  of  the  quarter  waveguide  at  the  KKKth  cross 
section.  The  fifth  page  of  RIVSWR  and  Fig.  7 
clearly  show  sides  1  through  4.  TANTi  is  nega¬ 
tive  (positive)  at  points  where  the  waveguide 
walls  slope  toward  (away  from)  the  z-axis. 

Read  in  from  EIGDAT.DAT 

This  file  saves  the  eigenvalue  and  coupling  coefficient 
data  for  the  taper.  It  can  be  used  to  rerun  the  code  for 
a  new  set  of  frequency  limits.  The  file  is  automatically 
read  when  a  "l”  is  typed  in  response  to  the  question, 
"EIGENVALUES  FOR  GUIDE  IN  EIGTAT .DAT?”  The  entries  of  a 
typical  row  are  defined  from  left  to  right  for  the  KKKth 


a.  COUPL(KKK) 


—  Solymar's  coupling  coefficient  (s^q)  [10]^ 

b.  USHS(KKK)  —  the  eigenvalue  squared  multiplied  by  the 

square  of  the  mesh  size 

c.  HH  (KKK)  —  mesh  size  (cm)  as  depicted  in  Fig.  5 

d.  FCUTOFF(KKK)  —  cutoff  frequency  normalized  by  l.E  +  9 

e.  ZLOC(KKK)  --  as  defined  for  RSIZ.DAT 

Output  Variables 

The  input  VSWR  and  scattering  matrix  are  output  by  the 
program.  The  following  is  a  list  of  key  variables  as  they  are 
defined  in  the  program. 

1 .  HZ  —  a  matrix  containing  of  the  current  cross 

section.  ""he  actual  field  points  are 
those  that  lie  on  and  interior  to  the 
solid  and  dashed  lines  of  Fig.  5. 

2.  FCUTOFF(KKK)  —  as  defined  for  EIGDAT.DAT 

3.  SI  1  ,S12, 

S21,S22  —  elements  of  the  transition  scattering 

matrix  at  a  specific  frequency 

4.  VSWR1  —  transition  input  VSWR  including  effects  of 

load  mismatch 

5.  VSWR2  —  transition  output  port  VSWR  including 

effects  of  source  mismatch 

6.  RATIO  —  cutoff  frequency  ratio  (ridged/unridged) 

at  the  current  cross  section  for  the 
dominant  TE  mode 
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Accuracy  Specifiers 

There  are  six  variables  that  can  be  used  to  determine  the 
accuracy  of  the  matrix-eigenvalue  problem  and  differential 
equation  solutions. 

1.  MAXITER  —  number  of  iterations  the  user  will  permit 

the  code  to  execute  in  an  attempt  to  obtain 
convergence  for  the  matrix-eigenvalue 
problem. 

2.  RESMAX  —  representative  of  the  maximum  allowed  error 

in  the  longitudinal  magnetic  field. 

3.  DELTAU  —  an  indicator  of  the  maximum  allowed  error  in 

the  cross-section  eigenvalue. 

4.  RTOL,ATOL  —  equivalent  to  RELERR  and  ABSERR  in  the 

differential  equation  solving  routine 
DESOLV.  They  are  representative  of  this 
routine's  accuracy. 

5.  KNT  —  this  variable  is  found  in  subroutine 

BLOCKS.  The  longest  dimension  (a  or  b)  is 
divided  into  KNT  equal  size  mesh  blocks  of 
length  and  width  h.  This  variable  is  found 
in  subroutine  BLOCKS.  KNT  is  incremented 
until  the  mesh  fits  the  waveguide  to  within 
RFIT.  In  the  current  code,  execution  is 
halted  at  KNT  =  1000  and  the  user  is  asked 
to  relax  the  tolerance  on  the  required  fit 
RFIT. 
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Error  Detection 


RIVSWR  responds  to  the  following  errors: 

1.  Required  Fit  Unobtainable  —  The  variable  IERR  is  returned 
from  subroutine  BLOCKS.  IERR  =  1  if  the  matrix  HZ  is  too 
small  to  obtain  a  block  size  capable  of  satisfying  the 
required  fit.  IERR  =  0  for  a  successful  fit.  If  IERR  = 
1,  the  user  is  notified  of  the  problem  and  is  prompted 
from  the  terminal  to  input  a  smaller  value  of  RFIT. 

2.  Frequency  Below  Cutoff  —  As  previously  mentioned, 
Solymar’s  normal  mode  analysis  is  invalid  at  and  below 
mode  cutoff.  The  current  propagation  frequency  is  com¬ 
pared  to  the  cutoff  frequency  of  each  cross  section.  If 
the  propagating  wave  condition  is  violated,  the  user  is 
notified  and  program  execution  stops  at  1000. 

A  flow  chart  for  RIVSWR  and  documented  listings  of  RIVSWR  and  DESOLV  are 
given  in  the  following  pages. 
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PROGRAM . - . . . RIVSWR 

AUTHOR - - - . - - - BRETT  BRAAT2 


PURPOSE:  TO  COMPUTE  THE  SCATTERING  MATRIX  AND  VSWR  FREQUENCY 

PROFILES  OF  AN  ARBITRARILY  TAPERED  DOUBLE  RIDGED-WAVEGUIDE. 

METHOD  i  THE  HELMHOLTZ  WAVE  EQUATION  IS  SOLVED  FOR  THE  DOMINANT 
MODE  MAGNETIC  FIELDS  Hz  AND  EIGENVALUES  AT  DISCRETE  AXIAL 
LOCATIONS  ALONG  THE  TAPER.  FINITE  DIFFERENCE  TECHNIQUES 
WERE  USED  TO  TRANSFORM  THE  HELMHOLTZ  EQUATION  INTO  A 
A  MATRIX-EIGENVALUE  PROBLEM.  THE  INVERSE  ITERATIVE  POWER 
METHOD  (SUCCESSIVE  OVER  RELAXATION)  IS  USED  TO  SOLVE  THIS 
PROBLEM.  THE  RESULTING  TE10  MODE  BOUNDARY  FIELDS  AND  THEIR 
DERIVATIVES  ARE  APPROXIMATED  AS  CONTINIOUS  FUNCTIONS  USING 
PIECEWISE  CONTINIOUS  CUBIC  SPLINES.  THE  SPLINES  ARE 
USED  TO  COMPUTE  SOLYMAR'S  COUPLING  COEFFICIENT  SCIOHIOI. 
THE  SELF  AND  IMPEDANCE  COUPLING  TERMS  ARE  CALCULATED 
USING  THE  EIGENVALUE  OF  THE  PARTICULAR  CROSS  SECTION. 

ALL  OF  THE  COUPLING  COEFFICIENTS  ARE  APPROXIMATED  AS 
CONTINIOUS  FUNCTIONS  OF  AXIAL  POSITION  USING  CUBIC 
SPLINES.  THE  COUPLED  SET  OF  DIFFERENTIAL  EQUATIONS  IS 
SOLVED  FOR  TWO  ORTHOGONAL  SETS  OF  MODE  AMPLITUDE  INITIAL 
CONDITIONS.  THE  SOLUTIONS  AT  Z  =  L  ARE  THEN  ALGEBRAICLY 
MANIPULATED  INTO  THE  FORM  OF  THE  RIDGED  TRANSITIONS 
SCATTERING  MATRIX.  THE  SCATTERING  MATRIX  AND  INPUT 
VSWR  ARE  WRITTEN  TO  SPARAM.DAT  AND  PVSWR.DAT,  RESPECTIVELY. 


SOLUTION  OF  HELMHOLTZ  EQUATION  WITH  NORMAL-GRADIENT  BOUNDARIES 
DOMINANT  TE  MODE  OF  RIDGE  WAVEGUIDE 


DIMENSION  HZ (1000, 1000) 

REAL  Y1 (1000) ,Y2( 1000) ,Y3( 1000) , Y4 ( 1000) , X 1 ( 1000) , X2 ( 1 000) , 

*X3 ( 1000)  , 

$X4 ( 1000) ,B( 1000) ,C (1000) ,D( 1000) ,TANT (4) , FREQUEN ( 200 ) .CQUPLdOl) , 
♦  USHS (101) t  HH  <101) .BETA (101) ,Z IMP (101) .FCUTOFdOl) .ZLOC(lOl) , 

$YY Y ( 4 ) , VSWR 1 (200) ,VSWR2(200) , T ANTI (101) ,TANT2(101) , TANT3 ( 101 ) , 

$T  ANT  4(101) 

REAL  EVAL , DERI V , U 

COMPLEX  All  .A12.A22.A21 , S 1 1 .S12.S21 , S22 , ZETA , GAMMA  1  .GAMMA2, 

IGAMAL , GAMAS 
INTEGER  I , N 
I NTE6ERM  MAX  I TR 
COMMON  /FVARS/  ZLOC , 

COMMON  NCROSS 
EXTERNAL  F 
MAX ITR  =  30000 
RESMAX  *  .001 
DELTAU  =  .0001 
SOL  =  2 . 997925E+ 1 0 
PI  =  3. 141593 

PERM  =  4 . *P I  *  1 . 0E-9  ‘MAGNETIC 
RELPER  =  1.0  ‘RELATIVE 


COUPL , Z I MP , BET  A 


‘ MAX  *  ITERATIONS  UNTIL  PROGRAM  QUITS 
IMAXIMUM  ACCEPTABLE  RESIDUE 
! MAX .  ACCEPT  RELATIVE  ERROR  IN  EIGENVALUE 
‘.SPEED  OF  LIGHT  CM/SEC 


PERM  I AB I L I T Y 
PERMITTIVITY 


OF 

OF 


FREE-SPACE (CM) 
FREE-SPACE 


0PEN(UNIT=16,FILE=  'SPAR AM. DAT ' , ST  AT US= ' UNKNOWN ' ) 
OPEN ( UN  I T  =  17, FILE* 'RSI Z. DAT ' , STATUS = ' UNKNOWN ' ) 
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OPEN (UNIT* 19 .FILE* 'PVSWR.DAT' ,  ST ATUS* ' UNKNOWN ' ) 
OPEN (UN I T*20, FILE*  EI6DAT.DAT' , ST ATUS* ‘ UNKNOWN ' ) 


C 

C  DATA  INPUT 
C 

WRITE(6,*) 'ENTER  REFLECTION  COEFFICIENT  OF  SOURCE' 

READ (A , *) GAMAS 

WRITE (6,*) 'ENTER  REFLECTION  COEFFICIENT  OF  LOAD' 

READ (6 , *) GAMAL 

WRITE (6,*) 1  ENTER  LOW  AND  HIGH  EDGES  OF  SWEEP  BAND  (GHz)' 

READ(A,*)FMIN,FMAX 

FREQUE  =  FMINM.OE  +  9 

WRI TE ! 6 , * ) '  ENTER  i  OF  FREQUENCY  STEPS' 

READ ( A , *) NFS 

DELF  =  < (FMAX-FNIN) /FLOAT (NFS) )  *1 . OE+9  ! STEP  SIZE  IN  FREQUENCY 

NFS  =  NFS  ♦  1  ! ONE  MORE  FREQ.  POINT  THAN  ft  FREQ.  STEPS 

WRITEI6,*) 'EIGENVALUES  FOR  TAPER  IN  E I GDAT . DAT?  TYPE  "1"  IF  SO' 

READ (6,*) I  YES 

IF ( IYES. EQ. 1 ) SO  TO  29 

WRITE (6 , 1 ) 

1  FORMAT ( IX , 'HOW  GOOD  SHOULD  THE  FIT  BE?  ( INCHES)  ' ) 

READ (6,*) RFIT 

READ (17,*) NCROSS 
DO  28  KKK  *  1, NCROSS 
KMM  =  KKK 

READ < 17 ,«)AAfBB,DD,SS, ZLOC (KKK) 

READ ( 17 , *) T  ANT  1 (KKK) , TANT2 (KKK)  ,TANT3(KKK) , T  ANT  4 (KKK ) 

ZLOC (KKK)  =  ZLOC (KKK) *2. 54 

2  CALL  BLOCKS (AA,BB,DD,SS,I BAR ,JBAR,IMAX,JMAX,WFIT,RFIT,H,IERR) 

IF ( IERR. EQ. 1 ) THEN 

WRITE(6,*) 'ACCURACY  ON  FIT  REQUIRES  A  MATRIX  LARGER  THAN  THIS 
$  PROGRAM  CAN  HANDLE:  TRY  A  LESS  ACCURATE  FIT  OR  CHANGE  THE 

*  HZ-MATRIX  AND  ALL  VECTORS  OF  THE  SAME  DIMENSION  (AT  THE 

$  BEGINNING  OF  THE  MAIN  PROGRAM)  TO  A  LARGER  SIZE.' 

WRITE  <6t 1) 

READ (6 , * ) RF I T 
GO  TO  2 
ELSE 
ENDIF 

WRITE (6,*) 'ERROR 
H  =  H*2. 54 
JMAX  =  JMAX+3 
IMAX  *  I  MAX  +2 
P  *  FLOAT (IMAX) 

Q  =  FLOAT ( JMAX ) 

C 

C  ALGORITHM  TO  DETERMINE 
C 

CINT  *  COS(PI/P) 

ALPHA  =  4.  /  (2.  -t-SQRT  (4.  -CINT*t2>  ) 

WRITE (6,*)  ACCELERATION  FACTOR  W  =  '.ALPHA 
B 1  =  FLOAT ( JMAX-3) 

Bll  =  B 1 *2 . 

A 1  =  FLOAT ( (IMAX  -  2) *2) 


OF  F I T=  ' , WF I T , '  (INCHES)  H  =  ',H,'  (INCHES)' 

!  CONVERT  MESH  SIZE  TO  CM 
I  CONVERT  ft  BLOCKS  TO  ft  NODES  ON  B/2 
!  CONVERT  ft  BLOCKS  TO  ft  NODES  ON  A/2 
!  ft  BLOCKS  IN  GUIDE  WIDTH 
!  ft  BLOCKS  IN  GUIDE  HEIGHT 

OPTIMUM  ACCELERATION  FACTOR 

+  CO=!PI/Q) 


BOA  =  Bll/Al 

UHSS  * ( P I / A1 ) **2  (RECTANGULAR  WAVEGUIDE  TE(IO) 

C 

C  NOTE:  FOR  THE  EACH  CROSS  SECTION  THE  APPROXIMATION 

C  USED  AS  A  START  I N6  VALUE  FOR  UHSQ  IS  THAT  OF  A  RETANGULAR 

C  WAVEGUIDE  HAVING  THE  SAME  SI2E. 

C 

C  ENTER  INITIAL  GUESS  FOR  FIELD  VALUES-  NOTE:  THE  INITIAL  GUESS 

C  FOR  THE  FIELD  VALUES  IS  ZERO  FOR  EACH  CROSS  SECTION. 

C 

DO  3  I  *  1 , 1  MAX 
DO  3  J  =  1 ,  J  MAX 
HZ ( I , J ) =0. 0 

3  CONTINUE 

DO  4  J  =  1 , JMAX-2 
HZ ( IMAX , J) =0. 

DO  4  I  —  1 1  IliAX  —  1 
HZ  (I ,J)=.5 

4  CONTINUE 

A2  =  FLOAT ( 2* J  BAR ) 

SOA  =  A  2  /  A 1 

B2  =  FLOAT (JMAX  -  IBAR  -  3) 

DINT  =  82*2. 

DOB  =  DINT / B 1 1 
C 

C  FIVE  ITERATIVE  CYCLES  FOLLOW 

C 

K0NVRG=0 

ITERAT=0 

5  ALFA=ALPHA/ <4. -UHSQ) 

DO  16  K=1 , 10 

C 

C  ITERATIVE  PASS  THROUGH  FIELD 

C 

BIGEST=0. 

IMA=  I  MAX  - 1 
J  =  1 

6  J=J*1 

IS  =  IMAX  -  JBAR 
JT  =  JMAX -  I  BAR- 1 
DO  7  1=2, IMA 

RESDL=HZ (1-1 , J)+HZ <1  +  1 , J) +  HZ ( I , J-l ) +HZ ( I , J  +  l ) ♦ (UHSQ-4. ) 
$  *HZ ( I , J ) 

IF (ABS (RESDL) .ST. BIGEST) B I GEST=ABS < RESDL ) 

HZ ( I ,J)=HZ  < I ,J) +ALFA+RESDL 

7  CONTINUE 

IF ( J - JT ) 6 , 8 , 9 

8  I  M  A  =  I  S 

HZ (ISM ,JT+1)=HZ  (IS-1  ,JT*1) 

9  IF  (J. LT. JMAX-1 ! GOTO  6 
C 

C  SET  HZ  VALUES  AT  EXTERNAL  FIELD  POINTS 
C 


HZ (1,1) =HZ (1,3) 

HZ ( 1 , JMAX ) =HZ ( I , JMAX-2 ) 

CONTINUE 

DO  11  I=IS*1 , IMAX-1 
HZ(I,JT+l)=HZ(I,JT-l) 

CONTINUE 
DO  12  J  =  1 , JHAX 
HZ ( 1 , J) =HZ (3 , J ) 

CONTINUE 

DO  13  J=JT+2 , JMAX 

HZ(IS+1,J)=HZ(I5-1,J) 

CONTINUE 

FIND  AVERA6E  HZ  NEAR  SYMMETRY  EDGE 

TOTAL*. 5* (HZ (IMAX-1 ,  2) *HZ ( IMAX-1 , JT) ) 

DO  14  J=3, JT-1 

TOTAL=TOTAL+HZ(IMAX-l,J) 

CONTINUE 
AVG=. 5*B2/T0TAL 

SCALE  FIELD  VALUES 

DO  15  I  =  1 , 1  MAX 
DO  15  J  =  1 , JMAX 

HZ ( I , J ) =AV6*HZ ( I ,  J  ) 

CONTINUE 
I TERAT  =  I TERAT  + 1 
CONTINUE 

END  OF  SINGLE  ITERATION 

CALCULATION  OF  RAYLEIGH  COEFFICIENT 

RCN=0 . 

RCD=0 . 

DO  17  J=2,JMAX-1 
DO  17  1=2, I  MAX  - 1 

FCTR  =  AREA(I , J, JT, IS, I  MAX, JMAX) 

RCN=RCN+FCTR*HZ < I , J ) * (HZ ( I  - 1 , J ) +HZ  < I  + 1 , J i +HZ ( I , J -  1 )  + 

1  HZ  (I ,  J  +  1J-4. *HZ (I  ,  J) > 

RCD=RCD+FCTR* (HZ ( I ,J)**2> 

CONTINUE 

RCN  =  RCN+.5*HZ (IS, JT+1 ) « (HZ ( IS-1 , JT+ 1 > -HZ ( IS+1 , JT-1  > ) 
RAYLGH=-RCN/RCD 

CHECK  FOR  CONVERGENCE 

IF  ( ABS ( ( R AYLGH-UHSQ ) /RAYLGH) . LT . DEL T AU ) KONVRG* 1 

IF  ( ITERAT.GT.MAXITR) GOTO  18 

UHSQ=R AYL6H 

IF (KONVRG.EQ.O)GOTO  5 

IF  (BIGEST.GT.RESMAX)GOTO  5 

FCUTOF  (K'KK)  =  (SORT  (UHSQ)  *SOL)  I  (  H» 2  .  « P I  *  1 .  OE  +  9 )  1  CUTOFF  FREQ  / 1 .  E<? 


7/  V  V  V  ’ 


.*  *.  >  .  ,  , 


WRITE  <6, *)  '  CUTOFF  FREQUENCY 


, FCUTOF ( KKK ) , '  GHz 


SCALE  THE  FIELD  VALUES  TO  SOLYMARS  NORMALIZATION 

DO  19  II  =  1 , 1  MAX 
DO  19  JJ  =  1 , JMAX 

HZIII.JJ)  =  HZ (II , JJ) / (SORT (4. *RCD*UHSQ) ) 

CONTINUE 

BELOW  IS  SHOWN  A  QUARTER  SECTION  OF  A  RIDGED  WAVEGUIDE.  THE 
LINES  OF  SYMMETRY  ARE  INDICATED  BY  DASHES  AND  THE  GUIDE  BOUNDARY 
IS  MADE  UP  OF  STARS.  THE  SECTION  BEING  LOOKED  AT  IS  THE  LOWER 
LEFT  ONE. 

FILL  UP  THE  X  AND  Y  VECTORS  FOR  SIDES  1-4  SHOWN  BELOW.  X  GIVES 
THE  BOUNDARY  LOCATION  AND  Y  IS  THE  FIELD  AMPLITUDE  AT  X.  THESE 
VECTORS  ARE  FIT  TO  A  CUBIC  SPLINE  AND  ARE  THEN  USED  IN  AN  INTE¬ 
GRATION  ROUTINE  THAT  COMPUTES  THE  BACKWARD  COUPLING  COEFF I E I  ENTS . 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 


SYMMETRY  PLANES  >>> 


<<<  SIDE  1  ( # I NT  =  J  MAX ) 


**************************  D/2 


SIDE  4  ( # I NT  =  J  BAR ) 


C  ♦  A  ! 

C  «  • I NTERVALS= I  MAX  SYMMETRY  PLANES  >>>  ! 

C  ♦  ! 

C  *  1 

C  * 

C  *  1 

C  *  **************************  D/2 

C  *  * 

C  *  <<<  SIDE  1  ( # I NT  =  J  MAX )  *  SIDE  4<#INT=JBAR> 

C  *  * 

C  *  SIDE  3  >>>  *<<<#  I NTERVALS= I  BAR 

C  *  * 

C  *  * 

C  B/2  ***********************************  A/2  -  S/2 

C 

C  SIDE  2 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

C  FILL  X  fc  Y  FOR  SIDE  1 


SIDE  3  >>>  *<<<#  I NTERVALS= I  BAR 


*•*»«»***«•*******«**«*****•«»***«  a/2  -  S/2 


DO  20  JP  =  2 , JMAX- 1 
J00=JP-1 

Xl(JOO)  =  (FLOAT (JOO-l )) *H 
Y 1 <J00)  =  HZ (2, JP) 

CONTINUE 

FILL  X  t,  Y  FOR  SIDE  2 

DO  21  IPP  =  2 , IS 
100  =  IPP-1 

X 2  ( 1 00 )  =  (FLOAT ( 100-1 ) >«H 
Y2 ( I  00 )  =  HZ ( IPP, JMAX-1 ) 
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CONTINUE 

FILL  X  l  V  FOR  SIDE  3 
KNT  =  1 

DO  22  JPP  =  JMAX- 1 , JT , - 1 
X3 (KNT)  =  < FLOAT  < KNT- 1 ) ) *H 
Y31KNT)  =  HZ  ( I  S , JPP ) 

KNT  =  KNTM 
CONTINUE 

FILL  X  t  Y  FOR  SIDE  4 
MCNT  »  1 

DO  23  IL  =  IS, I  MAX 

X  4 (HCNT)  =  (FLOAT (HCNT- 1 ) ) *H 
Y4 (HCNT)  =  HZ  < I L , J T ) 

MCNT  =  MCNT+1 
CONTINUE 

INTEGRATE  THE  SQUARED-DERI VI T IVES  OF  THE  FIELD  VALUES 
ON  EACH  SIDE  OF  THE  GUIDE  SHOWN  ABOVE 

COMPUTE  CONTRIBUTION  TO  SC10H103-  FROM  SIDE  1 

L  =  JMAX-2 

CALL  SPLINE(L,X1 ,Y1 ,B,C,D) 

GRAT1  =  0. 

DO  24  IR  *  1 , L  —  1 
BEG  *  XI (IR) 

END  *  XI (IR+1) 

CALL  GAUSS(L,BEG,END,X1  ,  Y 1  ,B,C,D,  A*'0) 

GRAT1  =  GRAT1  +  ANS 
CONTINUE 

GRAL1  *  GR AT  1  *TANT  1  K  .K) 

COMPUTE  CONTR'vUTION  TO  SC10K10I-  FROM  SIDE  2 

L  =  IMAX  -  JBAR  -  1 
CALL  SPLINE(L,X2,Y2,B,C,D) 

GRAT2  =  0. 

DO  25  IR  =  1 , L- 1 
BEG  =  X  2 ( I R ) 

END  =  X 2  <  I R -*■  1 ) 

CALL  GAUSS(L, BEG, END, X2,Y2.B,C,D, ANSI 
GRAT2  =  GRAT2  ♦  ANS 
CONTINUE 

GRAL2  =  6RAT2*TANT2  <K!-'K  > 

COMPUTE  CONTRIBUTION  TO  5  I  !  A ’  I  1 0  I  -  r'Z“  ZIZi  ' 

L  =  IBAR  ♦  1 

CALL  5PL I NE ( L , X  0 , Y  7 , P , C , D 1 


00  26  IR  -  i ,  I  BAR 
BE6  »  X  3 ( I R ) 

END  *  X3<  IR+1 ) 

CALL  GAUSS (L , BEG, END, X3 , Y3 , B , C , D , ANS) 

SRAT3  *  GRAT3  ♦  ANS 
CONTINUE 

BRAL3  *  6RAT3*TANT3 (KKK) 

COMPUTE  CONTRIBUTION  TO  StlOHlOl-  FROM  SIDE  4 
L  *  JBAR  ♦  1 

CALL  SPLINE(L,X4,Y4,B,C,D) 

BRAT4  *  0. 

DO  27  IR  *  1 ,  JBAR 
BE6  *  X 4 ( I R ) 

END  =  X4 ( IR+1 ) 

CALL  GAUSS (L, BEG , END ,  X4  ,  Y4 , B , C , D , ANS) 

GRAT4  =  GRAT4  +  ANS 
CONTINUE 

GRAL4  *  6RAT4*TANT4 (KKK) 

SUM  UP  ALL  FOUR  SIDES  AND  MULTIPLY  BY  4  TO  ACCOUNT  FOR  ALL  4 
QUADRANTS.  THEN  MULT  BY  -.5  TO  GET  SOLYMARS  COUPLING  COEFF. 

TINT  *  4. * (SRAL1+GRAL2+GRAL3+GRAL4) 

COUPL (KKK)  =  5*TINT  !  SOLYMARS  CDEFF 
USHS (KKK)  -  UHSQ 
HH  (KKK)  =  H 

SEND  CROSS  SECTION  COULPINB  COEFFICIENT  DATA  TO  EIGDAT.DAT. 

IF  THE  USER  SPECIFIES  THAT  EIGENVALUE  DATA  ALREADY  EXIST  THERE, 
THE  PROGRAM  SKIPS  ALL  OF  THE  ABOVE  CODE  AND  BEGINS  WORK  AT  LABEL 
29  BELOW. 

WRITE(20,*)C0UPL(KKK) ,USHS(KKK) ,HH (KKK) ,FCUTOF (KKK) , ZLOC (KKK) 
RATIO  =  PI/ (A1*SQRT (UHSQ) )  ! CUTOFF  FREQ-RATIO  WITH/WITHOUT 
CONTINUE 
GO  TO  31 

WRITE(6,*) 'ENTER  THE  NUMBER  OF  CROSS  SECTIONS  IN  EIGDAT.DAT' 

READ(6,*)NCR0SS 

DO  30  IMP  =  1 , NCROSS 

READ (20,*) COUPL (IMP) , USHS ( I MP) , HH ( IMP ) ,FCUTOF ( IMP) , ZLOC ( IMP) 
CONTINUE 
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c« ***»**•»*« #***#*#********#*****#»*#«*#**##*#*##*«***§* *» *#*#***###* **c 
c  c 

C  THIS  IS  A  SEPERATE  BLOCK  OF  THE  PROGRAM.  IN  THIS  SECTION,  THE  C 

C  TAPER  SCATTERING  MATRIX  AND  VSWR  ARE  COMPUTED  FROM  FMIN  TO  C 

C  FMAX  AT  INTERVALS  OF  DELF.  THE  ROUTINE  DESOLV  IS  USED  TQ  OBTAIN  C 

C  THE  AMPLITUDES  OF  A+  AND  A-  AT  Z=L  FOR  TWO  ORTHOGONAL  INITIAL  C 

C  CONDITION  VECTORS  DEFINED  AT  Z=0.  THESE  SOLUTIONS  ARE  THEN  C 

C  TRANSFORMED  INTO  THE  SCATTERING  MATRIX.  C 

C  C 

C  VSWR  AND  FREQUENCY  (NORMALIZED  BY  l.E+9  )  ARE  WRITTEN  TQ  THE  C 

C  DATA  FILE  PVSWR.DAT.  THE  SCATTERING  MATRIX  AND  FREQUENCY  ARE  C 

C  WRITTEN  TO  THE  FILE  SPARAM.DAT  IN  THE  ORDER  FREQUENCY , SI  1 , S22 ,  C 

C  S21.S22.  C 

C  C 

C****************»*************************************t***************C 

c 

31  DO  36  IFRQ  *  1 , NFS 

FREQUEN ( IFRQ)  =  FREQUE  ! IN  HERTZ  NOT  GHz 

WW  =  2. *PI*FREQUE 
DO  32  NBET  =  l.NCROSS 

CUTFRE  =  FCUTOF  <NBET)*1.0E+9 
IF (FREQUE. LE.CUTFRE) THEN 

WRITE (6,*)  'AT  THE  * ,NBET,  'th  CROSS  SECTION:  THE  CUTOFF 
$  FREQUENCY  IS  ' , CUTFRE,  '  Hz  ANT  T  'E  PROPAGATION  FREQUENCY 

*  IS  '.FREQUE,'  Hz.  START  OVER  WITH  HIGHER  PROPAGATION  FREQ' 

STOP  1000 


INITIALIZE  FORWARD  WAVE  TO  (1,0)  AT  Z=0 


UNITIALIZE  BACKWARD  WAVE  TO  (1,0)  AT  Z=0 


ELSE 
ENDIF 

BETA (NBET)  =  SQRT (RELPER# (WW/SOL ) **2-USHS (NBET) /HH (NBET > **2) 
Z IMP (NBET )  =  ALOG (WW*PERM/ BETA (NBET) ) 

CONTINUE 
NEN  *  4 
RTOL  *  l.OE-B 
ATOL  =  l.OE-B 
DO  35  I  *  1,2 
IF  <  I . EQ. 1 ) THEN 

YYY(l)  *  1.  !  INITIALIZE  FORWARD  WAVE  TO  (1,0 

YYY (2)  =  0. 

YYY (3)  *  0. 

YYY (4)  =  0. 

ELSE 

YYY ( 1 )  =  0.  UNITIALIZE  BACKWARD  WAVE  TO  (1,0 
YYY (2)  *  0. 

YYY (3)  *  1. 

YYY (4)  *  0. 

ENDIF 
TT  «  0.0 
TOUT  *  0.0 

DO  34  IKJ  «  1 , NCR0S5-1 
IFLAG  *  1 

TOUT  «  ZLOC  ( IKJM  ) 

CALL  DESOLV <F, NEN, YYY,TT, TOUT, RTOL, ATOL, IFLAG) 
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IF ( IFLA6. NE. 2) 60  TO  33 
IF ( I . EQ. 1 ) THEN 


All  > 
A21  • 
ELSE 
A12  > 
A22  * 
ENDIF 
CONTINUE 
CONTINUE 
ZETA  *  A1H 
Sll  *  -A21 
S22  *  All 


CHPLX ( YYY  <  1 )  , YYY (2) ) 
CHPLX (YYY (3) ,YYY (4) ) 


CNPLX < YYY  ( 1 ) , YYY ( 2) ) 
CHPLX (YYY (3) , YYY (4) ) 


RH0WE2) 


ZETA  *  A1 1  * A 2 2  -  A12*A21 

511  *  -A21/A22 

S22  *  A12/A22 

S21  *  ZETA/A22 

512  *  1/A22 

BAMHA1  =  SI  1*  ( (S12*S21*GAHAL) /  ( 1 . -S22*GAMAL) ) 
GAHMA2  *  S22+ ( (S12*S21 *GAMAS> / ( 1 . -SI  1*6 AM AS) ) 
RH0NE1  =  CABS ( GAMMA  1 ) 

RH0WE2  =  CABS (BAMHA2) 

VSMR1 (1FRQ)  =  (1.  +  RH0HE1 ) / ( 1 .  -  RH0NE1 ) 

VSHR2 ( IFRQ)  *  (1.  +  RH0WE2) / < 1 .  -  RH0WE2) 

THEFRE  =  FREQUE/ 1 . OE+9 

WRITE ( 16, *) THEFRE 

WRITE ( 16 , *) SI  1 , *  * , S 1 2 

WRITE ( 1 6 , *) S21 , *  * ,S22 

WRITE (19,*) THEFRE, *  * , VSWR1 ( IFRQ) 

FREQUE  =  FREQUE  +  DELF 
CONTINUE 
CLOSE ( 17) 

CLOSE (19) 

CLOSE (20) 

STOP 

END 


FUNCTION  SUBPR06RAM  TO  GENERATE  AREA  VALUES  AS  NEEDED 


INPUT  VARIABLES: 


I ,  J  -  LOCATION  OF  THE  POINT  TO  WHICH  AREA  VALUE  IS  FOUND 

IMAX.JMAX  -  SIZE  OF  THE  MATRIX  IN  I  AND  J,  RESPECTIVELY. 

JT , IS  -  POSITION  OF  THE  RIDGE  SIDE  AND  TOP 


OUTPUT:  AREA  -  INTEGRATION  AREA  ELEMENT  ASSOCIATED  WITH  THE  FIELD 
POINT  (I,J). 


FUNCTION  AREA(I,J,JT,IS,IMAX,JHAX) 
AREA=1 . 

IF  (1-2)30,1,2 
AREA*AREA/2. 

IF (J-JT) 12,8,3 
IF( I-IS) 5,4, 30 
AREA»AREA/2. 

IF(J-(JMAX-1))7,6,30 

AREA«AREA/2. 
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7  RETURN 

8  IF (I-IS)7,9, 10 

9  AREA*. 75 
RETURN 

10  AREA*AREA/2. 

U  IF  <1  —  1 MAX ) 7,6,30 

12  IF  t J-2) 30 ,10,11 
30  AREA*0. 

RETURN 

END 


SUBROUTINE  SPLINE  (N,  X,  Y,  B,  C,  D) 

INTEGER  N 

REAL  X (N) ,  Y (N) ,  B (N)  ,  C(N)  ,  D (N) 

C 

C  THE  COEFFICIENTS  B(I),  C(I),  AND  D(I),  1*1,2,... ,N  ARE  COHPUTED 
C  FOR  A  CUBIC  INTERPOLATING  SPLINE 
C 

C  S(X)  *  Y ( I )  +  B<I)«(X-X(in  +  C<I>*  CX-X  <1 ) )**2  +  D(I)*(X-X(I) )**3 
C 

C  FOR  X(I)  .LE.  X  .LE.  X(I+1) 

C 

C  INPUT.. 

C 

C  N  *  THE  NUMBER  OF  DATA  POINTS  OR  KNOTS  (N.GE.2) 

C  X  *  THE  ABSCISSAS  OF  THE  KNOTS  IN  STRICTLY  INCREASING  ORDER 

C  Y  *  THE  ORDINATES  OF  THE  KNOTS 

C 

C  OUTPUT.. 

C 

C  B,  C,  D  «  ARRAYS  OF  SPLINE  COEFFICIENTS  AS  DEFINED  ABOVE. 

C 

C  USING  P  TO  DENOTE  DIFFERENTIATION, 

C 

C  Y  ( I )  *  S  ( X  ( I ) ) 

C  B(I)  *  SP  (X  (I ) ) 

C  C(I)  *  SPP (X < I ) ) /2 

C  D ( I )  *  SPPP ( X ( I ) ) /6  (DERIVATIVE  FROM  THE  RIGHT) 

C 

C  THE  ACCOMPANYING  FUNCTION  SUBPROGRAM  SEVAL  CAN  BE  USED 
C  TO  EVALUATE  THE  SPLINE. 


INTEGER  NM1 ,  IB,  I 
REAL  T 

NM1  «  N-l 

IF  (  N  .LT.  2  )  RETURN 
IF  (  N  .LT.  3  )  GO  TO  50 

SET  UP  TRIDIAGONAL  SYSTEM 

B  *  DIAGONAL,  D  *  OFFDI A60NAL ,  C  *  RIGHT  HAND  SIDE. 


i 

ftt. 


u  o  u  u  u  u  u  u  u  u  u  u  u  u  u 


0 <  1 )  ■  X  (2)  -  X  ( 1 ) 

C (2)  «  <Y  (2)  -  Y < 1 ) ) /D ( 1 ) 

DO  10  I  *  2,  NM1 

D(l)  *  X(H-l)  -  X  ( I  > 

B ( I )  *  2- * (D < 1-1 )  ♦  D(I)) 

C  < 1  + 1 )  =  (Y ( 1+1 )  -  Y  < I ) ) /D ( I ) 

C ( I >  *  CCI+1)  -  C(I) 

10  CONTINUE 

END  CONDITIONS.  THIRD  DERIVATIVES  AT  X ( 1 )  AND  X(N) 
OBTAINED  FROM  DIVIDED  DIFFERENCES 

B ( 1 )  «  -Dtl) 

B(N)  *  -D(N-l) 

C ( 1 )  «  0. 

C(N)  *  0. 

IF  (  N  .EQ.  3  )  60  TO  15 
C  <  1 )  *  C (3) / (X (4) -X  (2) )  -  C (2) / < X (3) -X  ( 1 ) ) 

C <N)  *  C (N-l ) / (X (N) -X (N-2) )  -  C (N-2) / <X (N-l ) -X (N-3) > 

C < 1 )  *  C ( 1 ) «D ( 1 >  **2/ (X (4) -X ( 1 ) ) 

C(N)  *  -C(N)*D<N— 1)**2/(X(N>— X (N-3) > 

FORWARD  ELIMINATION 

15  DO  20  I  -  2,  N 

T  *  D(I-1)/B(I-1) 

B ( I )  =  B(I)  -  T*D  ( 1-1 ) 

C(I)  «=  C < I )  -  T *C < I  —  1 ) 

20  CONTINUE 

BACK  SUBSTITUTION 

C  CN)  «  C(N)/B(N> 

DO  30  IB  =  1,  NM1 
I  *  N-IB 

C(I)  =  (C(I)  -  D(I)*C(I+1) >/B(I> 

30  CONTINUE 

C ( I )  IS  NOW  THE  SIBMA(I)  OF  THE  TEXT 

COMPUTE  POLYNOMIAL  COEFFICIENTS 

B <N)  =  (Y  (N)  -  Y (NM1 ) ) /D (NH1 )  ♦  D (NM1 ) * (C (NM1 )  +  2.*C<N)> 
DO  40  I  *  1,  NM1 

B ( I )  *  ( Y ( 1  +  1 )  -  Y  < I ) > /D  < I )  -  D(I)»(C(I+1)  +  2.  *C (I ) ) 

D  ( I )  *  <cmi>  -  C  ( I ) )  /D  <  I ) 

C  ( I )  *  3 .  *  C  ( I ) 

40  CONTINUE 

C(N)  «  3 . * C ( N ) 

D(N)  *  D (N-l ) 

RETURN 
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SO  B ( 1 )  »  <Y(2)-Y<1) )/<X(2)-X(l) ) 
C<1)  ■  0. 

D(l)  «  0. 

B (2)  >  6(1) 

C(2)  -  0. 

D (2)  >  0. 

RETURN 

END 


SUBROUTINE  SEVAL (N,  U,  X,  Y,  B,  C,  D,  EVAL,  DERIV) 

INTE6ER  N 

REAL  U,  X (N) ,  Y(N),  BIN),  C(N),  D (N)  ,  EVAL,  DERIV 

THIS  SUBROUTINE  EVALUATES  THE  CUBIC  SPLINE  FUNCTION  AND  ITS 
FIRST  DERIVITIVE 

SEVAL  *  Y(I)  +  B(I)*(U-XU))  ♦  C  (I )  *  (U-X  (I ) )  **2  ♦  D  (I )  *  (U-X  ( I )  )  «*3 
DERIV  *  B(I>  +  2*C ( I ) *  <U-X (I ) )  ♦  3*D ( I > * (U-X (I > ) **2 

NHERE  X(I)  .LT.  U  .LT.  XU  +  l),  USIN6  HORNER'S  RULE 

IF  U  .LT.  XU)  THEN  I  -  1  IS  USED. 

IF  U  .BE.  X (N)  THEN  I  =  N  IS  USED. 

INPUT. . 

N  *  THE  NUMBER  OF  DATA  POINTS 

U  *  THE  ABSCISSA  AT  WHICH  THE  SPLINE  IS  TO  BE  EVALUATED 
X , Y  *  THE  ARRAYS  OF  DATA  ABSCISSAS  AND  ORDINATES 
B,C,D  *  ARRAYS  OF  SPLINE  COEFFICIENTS  COMPUTED  BY  SPLINE 

IF  U  IS  NOT  IN  THE  SAME  INTERVAL  AS  THE  PREVIOUS  CALL,  THEN  A 
BINARY  SEARCH  IS  PERFORMED  TO  DETERMINE  THE  PRDPER  INTERVAL. 

INTEGER  I,  J,  K 
REAL  DX 
DATA  1/1/ 

IF  (  I  . 6E.  N  )  I  «  1 

IF  (  U  .LT.  X(I)  )  60  TO  10 

IF  <  U  .LE.  XU  +  l)  )  BO  TO  30 

BINARY  SEARCH 

10  I  *  1 
J  *  N+l 

20  K  «  U+J)/2 

IF  (  U  .LT.  X (K)  )  J  B  K 

IF  (  U  .BE.  X (K)  )  I  *  K 

IF  (  J  .ST.  I+t  )  BO  TO  20 

EVALUATE  SPLINE 

30  DX  «  U  -  X(I) 
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EVAL  -  Y(I)  +  DX#(B(I)  ♦  DX*(C(I)  ♦  DX*Dd))) 

DERIV  *  B  ( I )  ♦  D  X  *  ( 2  *  C  ( I )  ♦  DX*3*D(I)) 

RETURN 

END 

C 

C  THIS  PROGRAM  COMPUTES  THE  INTEGRAL  OF  SOLYMAR'S  BACKWARD  C0UPLIN6 
C  COEFFICIENT  .  THE  METHOD  USED  IS  TO  PROVIDE  INTERVAL  END  POINTS 
C  (FROM, TO)  AND  THE  ENTIRE  ARRAY  OF  SPLINE  COEFFICIENTS  FOR  THE  FUNC- 
C  TION  BEING  INTEGRATED  (FUNCTION  IS  APPROXIMATED  BY  A  PEICEW1SE 
C  CONTINIOUS  CUBIC  SPLINE).  THE  SUBROUTINE  SEVAL  IS  USED  TO  EVALUATE 
C  THE  INTE6RAND  AT  THE  APPROPRIATE  GAUSSIAN  COORDINATES.  THE  VALUE  OF 
C  THE  INTEGRAL  BETWEEN  THE  END  POINTS  IS  -  ANSW 
C 

SUBROUTINE  GAUSS (L , FROM , TO , X , Y , B , C , D , ANSW) 

REAL  B(L) , C(L) ,D(L) , X (L) ,Y(L) , PS  I (4) , WEI (4) ,UV,VAL,DER 
DATA  PSI /-. 339901 ,.339901 ,-.861 136,. 061 136/ 

DATA  WEI/. 652145,. 652145,. 347855,. 347855/ 

C 

C  INPUT  VARIABLES: 

C 

C  FROM  -  LOWER  BOUND  OF  INTEGRAL  X(P) 

C  TO  UPPER  BOUND  OF  INTEGRAL  X(P+1) 

C  B,C,D  -  CUBIC  SPLINE  COEFFICIENTS  OBTAINED  FROM  ROUTINE-SPLINE 

C  L  LENGTH  OF  B,C,D  *  LENGTH  OF  X  OR  Y  IN  MAIN-LINE 

C  X , Y  -  DATA  POINTS  FOR  WHICH  SPLINE  COEFFICIENTS  WERE  FOUND 

C 

C  OUTPUT: 

C 

C  ANSW  -  VALUE  OF  THE  INTEGRAL 

C 

ANSW  *  0. 

DO  1  K  *  1,4 

UV  *  FROM  ♦  ( (TO-FROM) /2.)*(PSI(K)  ♦  1) 

CALL  SEVAL (L ,UV , X , Y ,B ,C , D , VAL , DER) 

FX  «  DER**2 

1  ANSW  *  ANSW*  WEI (K) *FX 

ANSW  =  ANSW« ( (TO-FROM) /2. ) 

RETURN 

END 

C 

C 

SUBROUTINE  BLOCKS (A, B,D,S, IBAR, JBAR, I  MAX , JMAX , WF I T , RF I T , H , I  ERR ) 

THIS  SUBROUTINE  COMPUTES  THE  THE  NUMBER  OF  BLOCKS  EACH  SIDE  OF  THE 
RIDGED  GUIDE  CROSS  SECTION  SHOULD  BE  BROKEN  UP  INTO.  AS  A  MINIMUM 
THE  LONGEST  SIDE  IS  DIVIDED  INTO  160  BLOCKS. 

THE  USER  SPECIFIES  A  REQUIRED  FIT,  IF  THE  WORST  FIT  ON 
ANY  OF  THE  GUIDE  SIDES  EXCEEDS  THIS,  THE  PROGRAM  WILL  ATTEMPT  TO 
INCREASE  THE  MESH  SIZE.  THIS  PROCESS  CONTINUES  UNTILL  THE  MESH 
IS  TOO  LARGE  (PRESENTLY  SET  TO  1000:  CAN  BE  MADE  LARGER)  OR  THE 
REQUIRED  FIT  IS  ACHIEVED. 


C 

C 


INPUT  VARIABLES: 


A,B|D|S 

< -  A . > 

««•**«»•  «••••**• 


******** 


••****« 


<  S  > 


OUTPUT  VARIABLES: 


IBAR,  JBAR 
I  MAX ,  JhAX 


REQUIRED  FIT  :  THIS  VARIABLE  IS  A 
LIMIT  ON  MAXIMUM  ERROR  OF  FIT  ON  ANY 
SUIDE  DIMENSION  (INCHES) 

STANDARD  RID6ED  WAVEGUIDE  DIMENSIONS 
(INCHES) 

A  -  WAVEGUIDE  WIDTH 
B  -  WAVEGUIDE  HEIGHT 
D  -  SPACE  BETWEEN  RIDGES 
S  -  RIDGE  WIDTH 


•  OF  BLOCKS  IN  RIDGE  HEIGHT  AND  HALF 

OF  *  IN  RIDGE  WIDTH 

HALF  THE  #  OF  BLOCKS  IN  GUIDE  WIDTH 

AND  HEIGHT,  RESPECTIVELY 

ACTUAL  FIT  OBTAINED,  WILL  ALWAYS  BE 

BETTER  THAN  REQUIRED  FIT (RFIT)  IF 

IERR  >  0 

1  IF  THE  MATRIX  HZ  IN  THE  MAIN 
PROGRAM  IS  TOO  SHALL  TO  FIT  TO  RFIT 
0  IF  RFIT  IS  MET 

THE  SIZE  OF  MESH  BLOCK  FOUND  TO 
SATISFY  RSIZ. 


CF  «  1.0/2. 

A2  «  A*CF 
B2  *  B*CF 
D2  *  D*CF 
S2  «  S*CF 
BMD  *  B2-D2 
WFIT  =  0.0 
KNT  *1Z>5 
KNT  *  KNT  +  1 
F  *  1. / (FLOAT (KNT) *2. ) 
IF  ( A2. LE. B2) THEN 
H  »  F«B 
ELSE 
fl  =  F*A 
ENDIF 

IF ( A2. LE. B2) THEN 
JMAX  «  KNT 
SIM  «  A2/H 
IMAX  «  NINT (SIM) 

ELSE 

I  MAX  «  KNT 
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SJH  *  B2/H 
JMAX  «  NINT(SJH) 

ENDIF 

SJB  =  S2/H 
J8AR  «  MINT (SJB) 

SIB  *  BMD/H 
IBAR  «  MINT (SIB) 

IF ( A2. LE. B2) THEN 

01  *  ABS (FLOAT ( IMAX ) *H  -  A2) 

ELSE 

01  *  ABS (FLOAT ( JMAX ) *H  -  B2) 

ENDIF 

02  *  ABS (FLOAT ( JBAR) *H  -  S2) 

03  =  ABS (FLOAT < IBAR) *H  -  BHD) 

NFIT  =  AMAX1 (01,02,03) 

IF (WFIT. 6T. RFIT) THEN 
IF(KNT.LT. 1000)60  TO  99 

IERR  *  1  ! NEED  MATRIX  LARGER  THAN  IN  MAIN  HZ  (1000, 1000) 
60  TO  999 
ELSE 

IERR  *  0 
ENDIF 
RETURN 
END 


EH 


$ 


SUBROUTINE  F  (T,  Y,YDOT) 

COMMON  /FVARS/  ZLOC,COUPL, ZIMP, BETA 
COMMON  NCROSS 

REAL  T, ZIMP ( 101 ) , BETA (101) ,COUPL( 101) ,EVAL,DERIV, 

$BBB ( 101 ) , CCC ( 101 ) , ODD ( 101 ) , ZLOC ( 101 ) , Y (4) ,YD0T(4) 

C* ##***#*##*#####*#**######*#*#***#*#**#*# ***** »«•«*«* ***************** 
C  SUBROUTINE  F  CONTAINS  THE  REAL  FORM  OF  SOLYMAR ' S  NORMAL  MODE  * 

C  EQUATIONS.  THESE  EQUATIONS  ARE  WRITTEN  IN  THE  FORMAT  REQUIRED  * 

C  BY  SHAMPINE'S  ORDINARY  DIFFERENTIAL  EQUATION  S0LVIN6  ROUTINE  * 

C  THERE  ARE  FOUR  EQUATIONS:  ONE  FOR  THE  REAL  AND  IMAGINARY  PART  OF  * 

C  THE  FORWARD  AND  BACKWARD  TEUOJ  NORMAL  MODES.  (SHAMPINES  ♦ 

C  ROUTINE  IS  CALLED:  DESOLV)  * 


3 


INPUT  VARIABLES: 


CURRENT  AXIAL  POSITION  WITHIN  THE  TRANSITION  (Z)  * 
A  VECTOR  CONTAINING  THE  AMPLITUDES  OF  THE  NORMAL  * 
MODES  Y(l)  *  A*  <  REAL ) ,  Y (2)  =  A+ ( IMAGINARY) ,  * 
Y (3)  =  A- (REAL) ,  AND  Y(4>  =  A- ( IMAGINARY) .  * 


OUTPUT  VARIABLES: 


C  YDOT  -  A  VECTOR  CONTAINING  dA/dz  AT  Z=T.  ITS  ELEMENTS  * 

C  ARE  ORDERED  LIKE  THOSE  OF  Y.  * 

C **»*•» ####*##**##♦#**############*#*##*** •*»****»*»****#*»»***»#»*#*#* 

c 

C  THE  PROPAGATION  CONSTANT  (B)  IS  A  FUNCTION  OF  Z  AND  IS  ONE  OF  THE 
C  COUPLING  COEFFICIENTS.  SUBROUTINE  SPLINE  FITS  A  CUBIC  SPLINE  TO 
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C  BETA  AND  GIVES  ITS  VALUE  (B)  AT  Z*T.  THIS  SAKE  PROCEEDURE  IS  USED 

C  TO  OBTAIN  SOLYMARS  COUPLING  COEFFICIENT  S.tlOJClO]  AND  dtln(IO]/dz 

C  AT  Z  «  T 
C 

C******************************* #****##*###*»**##*«#«*****»****»*»*##** 

c 

c 

CALL  SPLINE (NCROSS , ZLOC , BETA, BBB,CCC,DDD) 

CALL  SEVAL (NCROSS, T, ZLOC, BETA, BBB, CCC, DDD, EVAL, DERI  V) 

B  >  EVAL 

CALL  SPLINE (NCROSS, ZLOC, COUPL , BBB , CCC , DDD) 

CALL  SEVAL (NCROSS, T, ZLOC, COUPL, BBB, CCC, DDD, EVAL, DERI V) 

S  =»  EVAL 

CALL  SPLINE(NCRQSS,ZLQC,ZIHP,BBB,CCC,DDD) 

CALL  SEVAL (NCROSS, T, ZLOC ,Z IMP, BBB, CCC, DDD, EVAL.DERIV) 

Z  *  * . 5*DERI V 
D  *  S-Z 

YDOT(l)  =  B*Y (2)  +  D*Y (3) 

YDOT (2)  =  -B*Y ( 1 )  ♦  D*Y (4) 

YDOT (3)  =  D*Y <  1 )  -  B»Y (4) 

YDDT (4)  *  D*Y (2)  +  B*Y (3) 

RETURN 

END 
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C  DESOLV. FOR 

C  SHAMPINE  '  S  ODE  SOLVER 

C 

SUBROUTINE  DESOLV (F , NEQN , Y , T , TOUT , RELERR , ABSERR , I FLA6) 
LOGICAL  START, CRASH, STIFF 
DIMENSION  Y (NEQN) ,PSI (12) 

DIMENSION  YY(20) ,WT(20) ,PHI (20,16) ,P (20) ,YP<20) ,YPOUT (20) 

CC  COMMON  7CDE/  YY , WT , PHI ,P , YP , YPOUT , PSI 
EXTERNAL  F 

DATA  FQURU/2. 9802324e-B/ 

DATA  MAXNUM/500/ 

IF  (NEQN  .It.  1  .OR.  NEQN  .gt.  20)  BO  TO  10 
IF (T. EQ. TOUT )  BO  TO  10 

IF  (RELERR  .It.  0.0  .OR.  ABSERR  .It.  0.0)  GO  TO  10 
EPS  =  AMAX1 (RELERR, ABSERR) 

IF  (EPS  .le.  0.)  GO  TO  10 
IF  ( 1FLAS  .eq.  0.0)  GO  TO  10 
ISN  *  ISIGN(1 , 1  FLAG ) 

IFLAG  =  lABS(IFLAG) 

IF  (IFLAG  .eq.  1)  BO  TO  20 
IF  <T  .ne.  TOLD)  BO  TO  10 

IF  ( IFLAG  .ge.  2  .AND.  IFLAG  .le.  5)  BO  TO  20 
10  IFLAG  =  6 

RETURN 

20  DEL  *  TOUT-T 

ABSDEL  *  ABS (DEL) 

TEND  =  T+10. 0*DEL 

IF  (ISN  .It.  0)  TEND  =  TOUT 

NOSTEP  «  0 

KLE4  «  0 

STIFF  =  .FALSE. 

RELEPS  *  RELERR/EPS 
ABSEPS  =  ABSERR/EPS 
IF  (IFLAG  .eq.  1)  60  TO  30 
IF  (ISNOLD  .It.  0)  GO  TO  30 
IF  (DELS6N*DEL  .gt.  0.0)  GO  TO  50 
30  START  =  .TRUE. 

X  *  T 

DO  40  L  *  1 , NEQN 

40  YY (L)  *  Y (L ) 

DELSGN  =  SI6N ( 1 . 0 , DEL) 

H  *  SI6N(AMAX1 (ABS (TOUT-X) ,FDURU*ABS(X) ) ,TOUT-X) 

50  IF  (ABS(X-T)  .It.  ABSDEL)  GO  TO  60 

CALL  INTRP  (X,YY, TOUT, Y , YPOUT , NEQN , KOLD , PHI, PSI) 

IFLAG  =  2 

T  •  TOUT 
TOLD  «=  T 
ISNOLD  *  ISN 
RETURN 

60  IF  (ISN  .gt.  0  .OR.  ABS(TOUT-X)  .ge.  FOURU*ABS ( X ) )  60  TO  GO 
H  =  TOUT-X 
CALL  F ( X , YY , YP) 

DO  70  L  =  1 , NEQN 

Y (L)  *  YY (L) +H*YP (L) 
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ISN*5 


IFLA6  «  2 

T  «  TOUT 

TOLD  *  T 

ISNOLD  *  I5N 
RETURN 

IF  (NOSTEP  .It.  MAXNUM)  GO  TO  100 
IFLAS  =  ISN*4 
IF  (STIFF)  IFLAS  =  ISN*5 
DO  90  L  =  1 ,  NEQN 

Y  (L)  =  YY(L) 

T  *  X 

TOLD  =  T 

ISNOLD  =  1 

RETURN 

H  =  SI6N(AMIN1(ABS(H) ,ABS (TEND— X) ) ,H) 

DO  110  L  =  1 ,  NEQN 

HT  (L)  =  RELEPS*ABS(YY(L) )+ABSEPS 

CALL  STEP  <X,YY,F, NEQN, H,EPS,WT, START, HOLD, 
K,KOLD, CRASH, PHI  , P  ,  YP , PSI ) 

IF  (.NOT.  CRASH)  GO  TO  130 
IFLAS  =  ISN*3 
RELERR  *  EPS*RELEPS 
ABSERR  =  EPS*ABSEPS 
DO  120  L  =  1 , NEQN 

Y  (L )  *  YY  (L) 

T  =  X 

TOLD  =  T 
ISNOLD  =  1 

RETURN 

NOSTEP  *  NOSTEP+1 

KLE4  *  KLE4+1 

IF  (KOLD  .gt.  4)  KLE4  =  0 

IF  (KLE4  .ge.  50)  STIFF  =  .TRUE. 

GO  TO  50 
END 


SUBROUTINE  STEP  ( X , Y , F , NEQN , H , EPS , NT , START , HOLD , 

1  K, KOLD, CRASH, PHI ,P,YP, PSI) 

LOGICAL  START  , CRASH , PHASE  1 , NORND 

DIMENSION  Y (NEQN) , NT (NEQN) ,PHI (NEQN, 16) , P (NEQN) , YP (NEQN) , 
1  PSI (12) ,6STR(13) , TWO (13) 

DIMENSION  ALPHA (12) , BETA (12) , SIG ( 13) ,N ( 12) , V ( 12) , G ( 13) 
COMMON  /CSTEP/  GSTR , TWO , ALPHA , BETA , SI G , W , V , S 
EXTERNAL  F 

DATA  TWOU , FQURU/ 1 .4901162E-B,2.9B02324E-8/ 

DATA  TWO/2. ,4. ,8. , 16. ,32. ,64. , 1 28. ,256. ,512. , 

1  1024. ,2048. ,4096. ,8192./ 

DATA  GSTR  / . 5 , . 0833 , . 04 17 , . 0264 , . 0188 , . 0 t 43 , . 01 1 4 , 

1  .00936, .00789,. 00679,. 00592,. 00524,. 00468/ 
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jll&iAa  1 A -* &A.  S  hu, a 


v.<. 
Ia,‘  a  -t_X  -5  A 


DATA  8(1) ,6(2) , S IS ( 1 ) ✓ 1 . , . 5 , 1 . / 

CRASH  =  .TRUE. 

IF  (ABS(H)  .ge.  F0URU*A8S(X) )  GO  TO  5 
H  =  SIGN(FQURU*ABS (X) ,H) 

RETURN 

P5EPS  =  . 5*EPS 

ROUND  =  0. 

DO  10  L  =  1 , NEQN 

ROUND  =  ROUND+(Y(L)/WT(L) )««2 

ROUND  =  TNOU*SQRT (ROUND) 

IF  (P5EPS  .ge.  ROUND)  GO  TO  15 
EPS  *  2. * ROUND* ( 1 . +FOURU) 

RETURN 

CRASH  =  .FALSE. 

IF  (.NOT.  START)  SO  TO  99 
CALL  F  ( X , Y , YP) 

SUM  =  0. 

DO  20  L  =  1 , NEQN 

PHI (L , 1 )  =  YP (L) 

PHI (L , 2 )  *  0. 

SUM  =  SUM+(YP(L)/WT(L) )**2 
SUM  =  SORT (SUM) 

ABSH  =  ABS(H) 

IF  (EPS  .It.  If,.  *SUM*H*H)  ABSH  «  .  25*SQRT  (EPS/SUM) 

H  =  SIGN(AMAX1 <ABSH,FOURU*ABS(X) ) ,H) 

HOLD  =  0. 

K  =  1 

KOLD  =  0 

START  *  .FALSE. 

PHASE  1  =  .TRUE. 

NORND  =  .TRUE. 

IF  (P5EPS  .gt.  lOO.tROUND)  GO  TO  99 
NORND  =  .FALSE. 

DO  25  L  =  1 , NEQN 

PHI (L , 15)  =  0. 

IFAIL  =  0 

KPl  =  K  +  1 
KP2  =  K  +  2 
KM1  *  K-l 
KM2  =  K-2 

IF  (H  .ne.  HOLD)  NS  =  0 

NS  =  MINO (NS+ 1 , KOLD+ 1 ) 

NSP1  =  NS+1 

IF  (K  .It.  NS)  GO  TO  199 

BETA (NS)  =  1. 

REALMS  =  NS 
ALPHA (NS)  =  1. /REALMS 

TEMPI  =  H*REALNS 
SIS(NSPl)  =  1. 

IF  <K  .It.  NSP1 )  60  TO  110 

DO  105  I  =  NSP1,K 
IH1  =  1-1 

TEMP2  =  PSI <  IM1 ) 

PSI(IMl)  =  TEMPI 
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BETA  < I )  «  BETA ( IH1 ) *PSI <IMi)/TEMP2 

TEMPI  «  TEMP2+H 
ALPHA ( I )  *  H/TEMP1 

REALI  *  I 

S!6(I+1)  *  REALI* ALPHA ( 1 ) *SI 6  1 1 ) 

PSI  (K)  *  TEMPI 

IF  (NS  .gt.  1)  60  TO  120 

00  115  10  -  1  ,K 

TEMP3  *  IQ* ( IQ+1 ) 

V(IO)  *  1./TEMP3 
W (IQ)  -  V  C IB) 

60  TO  140 

IF  (K  .le.  KOLD)  60  TO  130 
TEMP4  *  K*KP1 

V  (K)  *  1./TEMP4 

NSM2  «  NS-2 

IF  (NSM2  .It.  1)  60  TO  130 

DO  125  J  =  1 , NSM2 

I  *  K-J 

V ( I )  «  V(I)-ALPHA(J*1)*V(I*1> 

LIMIT1  =  KP1-NS 
TEMP5  *  ALPHA (NS) 

DO  135  IQ  =  1 , LIMIT1 

V  ( IQ)  =»  V(IQ)-TEMP5*V(IQ*1) 

N ( IQ)  *  V< IQ) 

6 (NSP1 )  =  M(l) 

NSP2  *  NS+2 

IF  (KP1  .It.  NSP2)  60  TO  199 
DO  150  I  =  NSP2.KP1 
LIMIT2  *  KP2-I 
TEMP6  *  ALPHA ( 1-1 ) 

DO  145  IQ  «  1 , LIMIT2 

W ( IQ)  *  N(IQ)-TEMP6*N(IQ*1) 

I)  «  H  ( 1 ) 

CONTINUE 

IF  (K  .It.  NSP1)  60  TO  215 
DO  210  I  *  NSPl.K 
TEMPI  «  BETA  ( I ) 

DO  205  L  *  1 , NEQN 

PHI (L  ,  I )  =  TEMP1*PHI <L,I) 

CONTINUE 

DO  220  L  *  1 , NEQN 

PHI  (L ,KP2)  =  PHKLjKPl) 

PHI (L | KP 1 )  =  0.  !  TROUBLE? 

P (L)  «  0. 

DO  230  J  *  1 ,  K 

I  *  KP1-J 
IP1  «  1*1 

TEMP2  =  6(1) 

DO  225  L  =  1 , NEQN 

P (L)  ■  P (L) +TEMP2*PHI <L , I ) 

PHI (L  ,  I )  =  PHI  (L,I)+PHI  (L ,  IP  1  > 

CONTINUE 

IF  (NORND)  60  TO  240 


DO  23S  L  *  1 , NEON 

TAU  -  H*P (L) -PHI (L , 15) 

P (L)  *  Y (L) +TAU 

235  PHI ( L 1 16)  *  (P (L) -Y (L) ) “TAU 

60  TO  250 

240  DO  245  l  *  1 ,NEQN 

245  P(U  =  Y  (L)  +H*P (L) 

250  XOLD  «  X 

X  *  X+H 
ABSH  *  ABS(H) 

CALL  F  (X  ,P, YP) 

ERKM2  »  0. 

ERKM1  *  0. 

ERK  *  0. 

DO  265  L  *  I,  NEON 

TEMP3  *  l./HT(L) 

TEMP4  «  YP(L)-PHI(L,n 

IF  (KM2)  265,260,255 

255  ERKM2  *  ERKM2M  (PHI  <L , KM1  > +TEMP4)  *TEMP3)  **2 

260  ERKM1  *  ERKM 1  + ( <  PH I (L , K) +TEMP4 ) *TEMP3) **2 

265  ERK  =  ERK+  <TEMP4*TEMP3) **2 

IF  (KH2)  290,275,270 

270  ERKM2  *  ABSH*SI6 (KM1 ) *6STR (KM2) *SQRT (ERKM2) 

275  ERKM1  *  ABSH*SI6  (K)*6STR(KM1)*SQRT  ( ERKM  1 ) 

280  TEMPS  »  ABSH*SQRT  (ERK) 

ERR  *  TEMP5*(B(K)-B!KP1)) 

ERK  »  TEMP5*SI6(KP1)*6STR(K) 

KNEH  >  k 

IF  (KM2)  299,290,285 

285  IF  (AM AX  1 (ERKM1 , ERKM2)  .1*.  ERK)  KNEH  »  KM1 
BO  TO  299 

290  IF  (ERKH1  .le.  . 5*ERK)  KNEH  *  KM1 
299  IF  (ERR  .le.  EPS)  60  TO  400 
PHASE1  =  .FALSE. 

X  *  XOLD 

DO  310  I  =  l ,  K 

TEMPI  *  l . /BETA ( I ) 

1P1  »  1  +  1 

DO  305  L  *  i,NEQN 

305  PHI (L ,  I )  a  TEMPI* (PHI <L, I ) —PHI (L  ,  IP! ) ) 

310  CONTINUE 

IF  (K  .It.  2)  60  TO  320 

DO  315  I  =  2  ,K 

315  PSI(I-l)  =  PSI(I)-H 

320  IFAIL  »  IFAIL+1 

TEMP2  =  .5 

IF  (IFAIL— 3)  335,330,325 

325  IF  (P5EPS  .It.  . 25*ERK )  TEMP2  =  SORT (P5EPS/ERK) 
330  KNEH  a  i 

335  H  «  TEMP2*H 

K  >  KNEH 

IF  (ABS(H)  . ge.  FOURU*ABS(X>)  60  TO  340 
CRASH  a  .TRUE. 

H  *  S I  BN ( FOURU* ABS ( X ) ,H> 
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EPS  =  EPS+EPS 

RETURN 

SO  TO  100 

KOLD  «  K 

HOLD  «  H 

TEMPI  ■  H*6 (KP1 ) 

IF  (NORND)  60  TO  410 
DO  405  L  »  1 , NEQN 

RHO  *  TEMPI* (YP (L) -PHI (L,l))-PHI(L,16) 

Y  CL )  *  P (L) +  RHO 

PHI (L , 15)  *  (Y(L)-P(L) ) -RHO 

BO  TO  420 

DO  415  L  =  1 , NEQN 

Y (L )  *  P (L) +TEMP1* (YP(L)-PHI (Lfl) ) 

CALL  F  (X,Y,YP) 

DO  425  L  *  1 , NEQN 

PHI (L,KPl )  *  YP CL) -PHI  CL 1 1 > 

PHI (L , KP2)  =  PHI (L,KP1)-PHI (L,KP2) 

DO  435  I  *  1,K 

DO  430  L  =  1 , NEQN 

PHI (L , I)  *  PHI (L| I ) +PHI CL, KP 1 ) 

CONTINUE 
ERKP1  *  0. 

IF  (KNEW  .eq.  KM1  .OR.  K  .eq.  12)  PHASE1  * 

IF  (PHASED  GO  TO  450 

IF  (KNEW  .eq.  KM1 )  GO  TO  455 

IF  (KP1  .gt.  NS)  60  TO  460 

DO  440  L  -  1 , NEQN 

ERKP1  =  ERKP1 + (PHI <L , KP2) /WT (L) ) **2 

ERKP1  =  ABSH*GSTR (KP1 ) *SQRT (ERKP1 ) 

IF  (K  .gt.  1)  GO  TO  445 
IF  (ERKP1  .ge.  . 5*ERK)  GO  TO  460 
GO  TO  450 

IF  (ERKMl  .le.  AMIN1  (ERK.ERKPD)  GO  TO  455 
IF  (ERKP1  .ge.  ERK  .OR.  K  ,eq.  12)  GO  TO  460 
K  =  KP  1 
ERK  =  ERKP1 
GO  TO  460 
K  *  KM1 
ERK  *  ERKMl 
HNEW  *=  H+H 

IF  (PHASED  GO  TO  465 
if <p5eps. ge. erk*two (k+l ) )  go  to  465 
HNEW  «  H 

IF  (P5EPS  .ge.  ERK)  GO  TO  465 
TEMP2  *  K+l 

R  »  (P5EPS/ERK) **(1./TEHP2) 

HNEW  =  ABSH*AMAX1 (.5.AMIN1 C.9,R> ) 

HNEW  >  SIGN ( AMAX 1 (HNEW , FOURU* ABS ( X ) ) , H) 

H  =  HNEW 

RETURN 

END 


SUBROUTINE  INTRP  (X , Y , XOUT , YOUT , YPOUT ,NEQN, KOLD , PHI , PSI ) 
DIMENSION  Y (NEON) ,YOUT (NEON) , YPOUT (NEON) ,PHI (NEON, 16) 
DIMENSION  PSI (12) ,6(13) ,M(13) ,RH0(13> 

DATA  G(l) ,RHO(l)/l.  ,1./ 

HI  *  XOUT-X 
KI  =  KOID+1 
KIP  1  *  K I ♦ 1 

DO  5  I  »  1  *  KI 

TEMPI  =  I 

5  W ( I )  =  1. /TEMPI 

TERM  =  0. 

DO  15  J  =  2,KI 
JM1  =  J-l 
PSIJMt  =  PSI(JMl) 

GAMMA  «  (HI+TERM) /PSIJM1 
ETA  *  HI/PSIJH1 
LIMIT1  *  KIP1-J 
DO  10  I  «  l.LIMITl 
10  H(I)  *  BAMMA*N(I)-ETA*W(I+1) 

G(J)  =  H ( 1 > 

RHO ( 3 )  =  BAHMA*RHO ( JM1 ) 

15  TERM  *  PSIJM1 

DO  20  L  »  l.NEQN 
YPOUT (L)  *  0. 

20  YOUT (L)  *  0. 

DO  30  J  «  1 ,KI 

I  *  KIP1-J 
TEMP2  «  B ( I ) 

TEMP3  *  RHO ( I ) 

DO  25  L  »  1 , NEON 

YOUT (L)  *  Y0UT(L)+TEMP2*PHI(L,I) 

25  YPOUT (L)  *  YP0UT(L)+TEMP3*PHI (L,I> 

30  CONTINUE 

DO  35  L  *  l.NEQN 
35  YOUT  (L)  =  Y(D+HI*YOUT(L) 

RETURN 

END 


APPENDIX  C 


DATA  FILES  FOR  A  WR-90  TO  WRD-750  COSINE  IMPEDANCE  TAPER 

RSIZ.DAT 

This  file  contains  the  data  used  by  RIVSWR  to  compute  the  VSWR 
profile  of  the  cosine  impedance  transition  presented  in  Section  IV.  The 
first  line  contains  the  number  of  cross  sections  for  which  information 
is  given.  Every  two  lines  thereafter  contain  information  for  a  particu¬ 
lar  cross  section.  The  first  of  these  contains  a,  b,  d,  s  and  z,  which 
are  the  standard  dimensions  for  ridged  waveguide  and  axial  position, 
respectively.  The  second  line  contains  information  about  the  slopes  of 
the  waveguide  walls,  tan  01,  tan  02,  tan  03  and  tan  04.  Figure  7  in  the 
text  shows  how  they  correspond  to  the  waveguide  wall. 

EIGDAT.DAT 

This  file  contains  the  coupling  coefficient  data  which  were  saved 
by  RIVSWR. 

PVSWR.DAT 

The  first  and  second  columns  contain  frequency  (GHz)  and  input 
VSWR  data. 

SPARAM.DAT 

The  rows  of  this  data  file  contain  in  order, 

frequency  (GHz) 

Sll*  S12 

S21’  S22* 
etc . 
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/• 

L, 

s 


0.9000000 

0.4000000 

0.4000000 

0.1730000 

O.OOOOOOOE+OO 

-  • 

_  » 

-0.1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0. 1386580 

£ 

0.8979100 

0.3992100 

0.3966897 

0.  1730000 

1 . 0000001E-02 

E 

-0.1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0. 1465272 

0.8958200 

0.3984200 

0.3951645 

0.1730000 

2. 0000001E-02 

V 

r*J 

•  m 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0. 1538023 

0.8937300 

0.3976300 

0.3936478 

0. 1730000 

3. 0000001E-02 

:;'-j 

-0. 1045000 

-3. 9499998E-02 

0.0000000E+00 

-0. 1605045 

a 

0. 8916400 

0.3968400 

0.3886241 

0.  1730000 

4. 0000003E-02 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0. 1666548 

7? 

0.8895500 

0.3960500 

0.3834169 

0. 1730000 

5. 0000001E-02 

VJ 

-0. 1045000 

-3. 949999BE-02 

O.OOOOOOOE+OO 

-0. 1722735 

0. 8874600 

0.3952600 

0.3818693 

0.  1730000 

6. 0000002E-02 

•‘2 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0. 1773811 

a 

0.8853700 

0.3944700 

0.3780915 

0.  1730000 

7 . OOOOOOOE-02 

ri 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0. 1819977 

38 

0.8832800 

0.3936800 

0.3724256 

0. 1730000 

8.0000006E-02 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0. 1861429 

b  - 

0.8B11900 

0.3928900 

0.3710971 

0.  1730000 

9. 0000004E-02 

\d 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0. 1898363 

0.8791000 

0.3921000 

0.3654709 

0. 1730000 

0. 1000000 

a 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0. 1930969 

w 

0.8770100 

0.3913100 

0.3623646 

0.  1730000 

0.1100000 

V\ 

> 

-0.1045000 

-3.9499998E-02 

O.OOOOOOOE+OO 

-0. 1959438 

V 
■  * 

0.8749200 

0.3905200 

0.3578277 

0. 1730000 

0. 1200000 

«r 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0. 1983955 

£ 

0.8728300 

0.3897300 

0.3541512 

0. 1730000 

0. 1300000 

8 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0.2004704 

0.8707400 

0.3889400 

0.3495793 

0.1730000 

0. 1400000 

>; 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0.2021862 

0. 8686500 

0.3881500 

0.3463731 

0. 1730000 

0.1500000 

A 

" 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0.2035608 

*2 

0.8665600 

0.3873600 

0.3408964 

0. 1730000 

0. 1600000 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0.2046114 

ty 

. * 

0.8644700 

0.3865700 

0.3383701 

0. 1730000 

0. 1700000 

$ 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0.2053549 

0.8623800 

0.3857800 

0.3329767 

0. 1730000 

0. 1800000 

%■ 

-0. 1045000 

-3. 949999BE-02 

O.OOOOOOOE+OO 

-0.2058081 

s." 

It* 

0.8602900 

0.3849900 

0.3291577 

0. 1730000 

0.1900000 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0.2059870 

-W 

0.8582000 

0.3842000 

0.3261495 

0.1730000 

0.2000000 

* 

■ 

-0. 1045000 

-3. 949999BE-02 

O.OOOOOOOE+OO 

-0.2059076 

v. 

0.8561100 

0.3834100 

0.3206090 

0. 1730000 

0.2100000 

$ 

-0. 1045000 

-3. 9499998E-02 

O.OOOOOOOE+OO 

-0.2055855 

0.8540200 
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2.362200 

3. 7761207E-03 

4. 7530641E-05 

5. 3825048E-03 

6.111437 

2.387600 

3. 7853597E-03 

4. 7530641E-05 

5. 3665149E-03 

6.  129646 

2.413000 

3. 7965600E-03 

4. 7530641E-05 

5.3505250E-03 

6. 147964 

2.438400 

3. 8078506E-03 

4. 7530641E-05 

5. 3345351E-03 

6.166392 

2.463800 

3. 8192752E-03 

4. 7530641E-05 

5 . 3 1 85457E-03 

6.  184930 

2.489200 

3. 8307908E-03 

4. 7530641E-05 

5. 3025563E-03 

6.2035B0 

2.514600 

3. 6964428E-03 

4. 7315221E-05 

5. 2B65660E-03 

6.208228 

2.540000 

-  110  - 


8 


PVSMR.DAT 

8.400000 

1.386806 

8.495999 

1.360496 

8.592000 

1.335621 

8.688000 

1.312075 

8.784000 

1.289767 

8.880000 

1.268615 

8.976000 

1.24854B 

9.072000 

1.229506 

9. 168000 

1.211434 

9.264000 

1.194281 

9.360000 

1.177998 

9.456000 

1.162550 

9.552000 

1.147899 

9.648000 

1.134008 

9.744000 

1.120849 

9.840000 

1.108394 

9.936000 

1.096617 

10.03200 

1.085500 

10.12800 

1.075025 

10.22400 

1.065181 

10.32000 

1.055962 

10.41600 

1.047380 

10.51200 

1.039461 

10.60800 

1.032271 

10.70400 

1.025941 

10.80000 

1.020730 

10.89600 

1.017095 

10.99200 

1.015596 

11.08800 

1.016362 

11.18400 

1.018780 

11.28000 

1.022060 

11.37600 

1.025675 

11.47200 

1.029344 

11.56800 

1.032919 

11.66400 

1.036318 

11.76000 

1.039505 

11.85600 

1.042452 

11.95200 

1.045149 

12.04800 

1.047591 

12. 14400 

1.049772 

12.24000 

1.051700 

12.33600 

1.053375 

12.43200 

1.054802 

12.52800 

1.055988 

12.62400 

1.056941 

12.72000 

1.057667 

12.81600 

1.058175 

12.91200 

1.058472 

13.00800 

1.058569 

13. 10400 

1.058473 

13.20000 

1.058195 

13.29600 

1.057744 

.-v  v 

*  "  V  -V 

13.39200 

1.057128 

13.48800 

1.056356 

13.58400 

1.055440 

13.68000 

1.054389 

13.77600 

1.053209 

13.87200 

1.051913 

13.96800 

1.050509 

14.06400 

1.049005 

*  14.16000 

1.047411 

!  14.25600 

1.045736 

1  14.35200 

1.043987 

1  14.44800 

1.042176 

J  14.54400 

1.040308 

^  14.64000 

1.038393 

14.73600 

1.036438 

14.83200 

1.034452 

I  14.92800 

1.032442 

1  15.02400 

1.030416 

15.12000 

1.028383 

15.21600 

1.026349 

1  15.31200 

1.024323 

15.40800 

1.022313 

15.50400 

1.020329 

15.60000 

1.018381 

l  15.69600 

1.016481 

15.79200 

1.014645 

15.88800 

1.012891 

15.98400 

1.011249 

1  16.08000 

1.009754 

16. 17600 

1.00B466 

16.27200 

1.007464 

16.36800 

1.006B38 

16.46400 

1.006659 

16.56000 

1.006924 

1  16.65600 

1.007553 

1  16.75200 

1.008432 

16.84800 

1.009464 

16.94400 

1.010577 

17.04000 

1.011725 

17.13600 

1.012876 

1  17.23200 

1.014010 

I  17.32800 

1.015115 

17.42400 

1.016177 

17.52000 

1.017191 

17.61600 

1.018152 

17.71200 

1.019056 

17.80800 

1.019900 

;  17.90400 

1.020680 

18.00000 

i 

1.021397 

SPARAM.DAT 
B. 400000 

(0.1117673,0.1173524)  (-0.9621535,-0.2190810) 

(-0.9621542,-0.2190819)  (-0.  1515571 ,5. 739409  IE-02) 

8.495999 

(0.1125035,0. 1032788)  (-0.9777583,-0.1437519) 

(-0.9777579,-0.1437549)  (-0. 1374689,6. 652B171E-02) 

8.592000 

(0. 1 12001 0,9. 0025298E-02)  (-0. 9872565 , -6. 8376124E-02) 

(-0. 9872574,-6. 8378299E-02)  (-0.  1233421,7. 372605  IE-02) 

8.688000 

(0. 1104124.7.7637993E-02)  ( -0. 9908267 , 6. 6189 100E-03) 

(-0. 9908268,6. 6208006E-03)  (-0. 1093659 , 7 . 91 06092E-02) 

8.784000 

(0. 1078841,6.614B028E-02)  (-0. 988661 3 , 8. 0B318B5E-02 ) 

( -0.9886622, 8. 0832034E-02)  (-9. 5707238E-02 , 8 . 2793549E-02) 

8.880000 

(0. 1045513,5. 5576466E-02)  (-0.9809687,0. 1538844) 

(-0.9809697,0. 1538B30)  ( -8. 251 1 701 E-02 , 8. 49206B9E-02) 

8.976000 

(0. 1005426, 4. 5930579E-02)  (-0.9679709,0.2254196) 

(-0.9679708,0.2254197)  (-6. 9906831 E-02, 8. 5624255E-02) 

9.072000 

(9.5980011 E-02, 3. 7209637E-02)  (-0.9498990,0.2951175) 

(-0.9498993,0.2951186)  (-5. 7996B43E-02 , 8. 50476 1BE-02) 

9. 168000 

(9.0976208E-02,2. 9402930E-02)  (-0.9270008,0.3626694) 

(-0.9270009,0.3626684)  (-4. 6B71435E-02 ,8. 3332255E-02) 

9.264000 

(8.5635200E-02,2.2491097E-02)  (-0.8995230,0.4278067) 

(-0.8995221,0.4278093)  ( -3. 6594670E-02 , 8. 0623366E-02) 

9.360000 

<8. 0052875E-02,1. 644976  IE-02)  (-0. B677302 , 0. 49027 1  1  ) 

(-0.8677298,0.4902718)  (-2. 7219538E-02 , 7. 7059306E-02) 

9.456000 

(7.4319802E-02,!. 1247266E-02)  (-0.8318903,0.5498251) 

(-0.8318903,0.5498261)  (-1 . B7B2 1 52E-02 , 7 . 27B1B91E-02) 

9.552000 

(6. 85 1 6247E-02 , 6, 8462300E-03)  (-0. 7922795,0.6062606) 

(-0.7922778,0.6062622)  ( - 1  .  1301 333E-02 , 6. 7923509E-02) 

9.648000 

(6. 2714346E-02 ,3. 2056409E-03)  (-0.7491630,0.6594013) 

(-0. 7491652,0.6593999)  ( -4. 78 1 4809E-03 , 6. 261 4240E-02 ) 

9.744000 

(5.6980848E-02,2. 81  1909  IE -04)  (-0.7028353,0.  7090666) 

(-0.7028392,0.7090626)  ( 7 . 84 1 2448E-04 , 5 . 6975521 E-02) 

9.840000 

(5.  1372845E-02.-  1 . 9736623E-03 )  (-0.6535717,0.7551  165) 

(-0.6535684,0.7551183)  (5. 4 1 43695E-03 , 5 . 1 124550E-02) 

9.936000 

<4.5940988E-02,-3.6091241E-03)  (-0.6016572,0.7974241) 

(-0.6016585,0. 797422B)  ( 9. 1 3956 1 l E-03 , 4 . 5 1 67 1 22E-02 ) 

10.03200 

<4.0729750E-02,-4. 6761 045E-03 )  (-0. 5473729,0. 8358838) 
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(-0.5473739,0.8358832)  ( 1 . 199577BE-02 , 3. 9202962E-02) 

10. 12800 

(3. 57767 4 2E -02 , -5. 2250270E-03)  (-0.4910094,0.8704033) 

(-0.4910096,0.8704034)  ( 1 . 4031467E-02 , 3. 332261 7E-02) 

10.22400 

(3.11 1 1963E-02 , -5. 3094020E-03)  (-0.4328377,0.9009192) 

(-0.4328385,0.9009186)  ( 1 . 529751 6E-02 ,2. 760675SE-02) 

10.32000 

(2. 6760191E-02 , — 4. 979449BE-03)  (-0.3731494,0.9273719) 

(-0.3731485,0.9273723)  ( 1 . 585381  IE-02 , 2 . 21 26 18 1 E-02) 

10.41600 

(2.2741562E-02.-4.2865803E-03)  (-0.3122055,0.9497324) 

(-0.3122081,0.9497321)  < 1 . 5762232E-02, 1 . 6943250E-02) 

10.51200 

( 1 . 9068B22E-02 , -3.2807277E-03)  (-0.2503003,0.9679750) 

(-0.2502979,0.9679755)  ( 1 . 5088214E-02, 1 . 21 13265E-02) 


10.60800 

(1. 5751589E-02,— 2. 01041B3E-03)  (-0.  1876745,0.9821028) 

(-0. 1876734,0.9821035)  ( 1 . 39002 14E-02 , 7 . 6764380E-03) 

10.70400 

( 1 . 27937 60E -02 ,-5.2206113E-04)  (-0. 1246102,0.9921234) 

(-0. 1246101,0.9921231)  < 1 . 22671 12E-02 , 3. 6693804E-03) 

10.80000 

(1. 0194955E— 02, 1 . 1395689E-03)  (-6. 1 360449E-02 , 0. 9980632) 

(-6.  1358325E-02, 0.9980632)  ( 1 , 0258 1 46E-02 , 1 . 178701 6E-04 ) 

10.89600 

(7.9516694E-03.2.9320165E-03)  ( 1 . B257959E-03 , 0. 9999622 ) 

(1.8280762E-03, 0.9999624)  (7 . 94 1 0020E-03 , -2. 9608773E-03 ) 

10.99200 

(6.0564550E-03,4.B15B932E-03)  (6. 4708143E-02 , 0. 9978740) 

(6. 4709552E-02, 0.9978738)  <5. 3B38054E-03,-5. 5576507E-03) 

11.08800 

(4.4990922E-03.6.7532999E-03)  (0.1270497,0.9918628) 

(0. 1270450,0.9918640)  <2. 6514668E-03 , -7 . 6690544E-03) 

1 1 . i 9400 

<3.2662388E-03,8.7106302E-03)  (0. 1 8860B2 , 0. 9820082) 

<0. 1886066,0.9820088)  < - 1 . 9339587E-04 , -9 . 3008662E-03 ) 

11.28000 

(2.3419673E-03, 1.0655509E-02)  (0.2491723,0.9683977) 

(0.2491721,0.9683977)  ( -3. 09 1 9253E-03 ,- 1 . 0462451 E-02) 

11.37600 

(1. 708402 1 E-03, 1 . 2559048E-02 )  (0.3085276,0.9511308) 

(0.3085228,0.9511328)  ( -5. 9892759E-03 , - l . 1 1 70791 E-02) 

11.47200 

(1.3473446E-03.1.4396776E-02)  (0.3664681,0.9303185) 

(0.3664690,0.9303181)  ( -8 . 8333 1 06E-03 , - 1  .  1 448 1 03E-02 ) 

11.56800 

( 1. 237 1755E-03, 1.61 45406E-02)  (0.4228041,0.9060763) 

(0.4228050,0.9060761)  (-1.  1 57B434E-02 , - 1  .  1 31 9584E-02) 

11.66400 

( 1.3571 01 0E-03, 1 . 7  783500E-02)  (0.4773310,0.8785419) 

(0.4773280,0.8785439)  < - 1 . 4 1 8 1 690E -02  ,  -  1 . 08 1 5879E-02 ) 

11.76000 

(1. 6832743E-03, 1 . 9296775E-02)  (0.5298963,0.8478417) 


(0.5298942,0.8478423)  <- 1 . 6607683E-02 , -9. 9692577E-03) 

11.85600 

(2. 19427568-03, 2. 0668905E-02)  (0.5803131,0.8141282) 

(0.5803130,0.8141285)  (- 1 . BB22B58E-02 , -8. 81 60457E-03) 

11.95200 

(2. 8645562E-03 ,2. 1889631E-02)  (0.6284368,0.7775474) 

(0.6284371,0.7775475)  <-2.0801B42E-02,-7.3920032E-03) 

12.04800 

<3. 67323988-03,2. 2950221 E-02)  (0.6741149,0.7382606) 

(0.6741119,0.7382634)  (-2 . 2522807E-02 , -5. 73B8390E-03) 

12.14400 

(4.59593128-03,2. 384294 0E -02)  (0.7171994,0.6964446) 

(0.7171991,0.6964450)  (-2. 3967544E-02,-3. B937540E-03) 

12.24000 

(5. 61023928-03,2. 4566106E-02)  (0. 7575B62,0. 6522481) 

(0. 7575880,0.6522468)  (-2.5127161E-02.-1 . 8974184E-03) 

12.33600 

(6.6939294E-03,2.51 16965E-02)  (0.7951410,0.6058667) 

(0.7951396,0.6058693)  (-2. 5992554E-02.2. 1 105599E-04) 

12.43200 

(7.8269243E-03,2.5495838E-02)  (0.8297742,0.5574611) 

(0.8297738,0.5574615)  (-2. 6562475E-02 , 2. 3918946E-03) 

12.52B00 

<8.9B80B80E-03,2.5705658E-02)  (0.8613789,0.5072331) 

(0.8613781,0.5072348)  ( -2. 683878  IE-02 , 4 . 6084910E-03 ) 

12.62400 

( 1 . 0158765E  —  02 , 2 . 5750922E  —  02)  (0.8898757,0.4553622) 

(0.889B74B, 0.4553646)  ( -2. 6828013E-02 , 6. B242606E-03) 

12.72000 

(1. 1321092E-02.2. 5637235E-02)  (0.9151924,0.4020407) 

(0.9151925,0.4020411)  (-2. 6539642E-02 , 9 . 0051256E-03) 

12.81600 

( 1 . 2 457 46  IE -02 ,2.5372019E-02)  (0.9372662,0.3474669) 

(0.9372661,0.3474671)  (-2. 5985766 E -02,1. 11 19625E-02) 

12.91200 

(1.3553559E-02,2. 4963653E-02)  (0. 95604B9 , 0. 2918290) 

(0.9560509,0.2918226)  < -2. 5184194E-02 , 1 . 3 139052E-02) 

13.00800 

( 1.4595901 E-02, 2. 44221 15E-02)  (0.9715022,0.2353172) 

(0.9715041,0.2353090)  ( -2 . 4153737E-02 , 1 . 5035833E-02) 

13.10400 

( 1. 557 1354E-02, 2. 3757981 E-02)  (0.9835988,0.  1781195) 

(0.9835985,0.1781229)  ( -2 . 29 1 39 10E-02 , 1 . 678BB 12E-02) 

13.20000 

<1.646B560E-02,2.29B3B23E-02)  (0.9923156,0. 1204580) 

(0.9923162,0. 1204545)  ( -2. 1 4B9425E-02 , 1 . 8376049E-02) 

13.29600 

(1.727971 1E-02,2.2110673E-02)  (0. 9976493,6. 251 6995E-02) 

(0. 9976496,6. 25 13031E-02)  (-1 . 9904645E-02, 1 . 9780476E-02) 

13.39200 

(1.7995946E-02.2. 1150874E-02)  (0 . 9996044 , 4 . 4866577E-03 ) 

<0. 9996043,4. 4858246E-03)  (- 1 . 8 1 84863E-02 , 2. 0988667E-02 ) 

13. 48800 

<  1 . 8610945E-02 , 2. 01  17678E  —  02)  (0. 9981942,-5.  3450249E-02) 
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( -1 . 6356045E-02, 2.  1 990042E-02) 


(0.9981948, 
13.58400 
(1.91 1975BE 
<0. 9934462, 
13.68000 
( 1 . 9520009E 
(0.9853922, 
13.77600 
( 1 . 9808996E 
(0.9740784, 
13.87200 
( 1 . 9985864E 
(0.9595637, 
13.96800 
(2. 0051392E 
(0.9419074, 
14.06400 
(2. 0007469E 
(0.9211857, 
14. 16000 
( 1 . 9856857E 
(0.8974825, 
14.25600 
( 1 . 9603299E 
(0.8708873, 
14.35200 
( 1 . 9251  BABE 
(0.8414994, 
14.44800 
( 1 . 8808600E 
(0.8094350, 
14.54400 
( 1 . 8279135E 
<0.7747915, 
14.64000 
(1. 7670654E 
(0.7377083, 
14.73600 
( 1 . 699087BE 
(0. 6983079, 
14.83200 
( 1 . 6247939E 
(0.6567296, 
14.92800 
( 1 . 5450008E 
<0.6131012, 
15.02400 
( 1 . 4605707E 
(0.5675898, 
15. 12000 
( 1 . 372361  IE 
(0.5203249, 
15.21600 
(1.281 2630E 


-5. 3449791E-02) 

-02,1. 9024707E-02)  (0.9934451,-0.  1110824) 

-0.1110792)  (-1 . 4445409E-02,2. 2777800E-02) 

-02,1.7BB4593E-02)  (0.9853925,-0. 1682300) 

-0.  1682298)  (-1 . 2480588E-02, 2. 3347927E-02) 

-02 , 1 . 6709087E-02)  (0.9740784,-0.2247203) 

-0.2247225)  (-1 . 048675  IE-02 , 2. 3698302E-02) 

-02.1.5512B03E-02)  (0.9595641,-0.2003511) 

-0.2803516)  <-B. 490458  IE-03, 2. 3832725E-02) 

-02, 1.4306905E-02)  (0.9419075,-0.3349680) 

-0.3349677)  (-6. 5154228E-03 ,2. 3754729E-02) 

-02, 1 . 3103309E-02)  (0.9211859,-0.3883863) 

-0.3883878)  (-4. 5866026E-03 , 2. 34723 17E-02) 

-02,1.  1913785E-02)  (0.8974022,-0.4404422) 

-0.4404424)  (-2. 7246689E-03 , 2. 2995798E-02) 

-02,1.0748558E-02)  (0.8708872,-0.4909744) 

-0.4909736)  (-9. 5143204E-04,2. 23365B2E-02) 

-02,9.6174562E-03)  (0. B4 14994, -0.5398288) 

-0.5398285)  (7. 1568845E-04,2. 1508692E-02) 

-02, B. 5303327E-03)  (0. 8094351 ,-0. 5B6B465) 

-0.5868466)  (2. 2593676E-03 , 2. 0528750E-02) 

-02,7.4948603E-03)  (0.7747916,-0.6319085) 

-0.6319090)  (3. 666524  IE- 03, 1 . 9412B34E-02) 

-02, 6. 519731  IE  — 03)  (0.  7377092,-0.6748545) 

-0.6748569)  ( 4. 9243635E-03 , 1 . 81 80005E-02) 

-02,5. 6103249E-03)  (0. 6983085,-0.7155725) 

-0.7155737)  (6.0232584E-03, 1 . 6849030E-02) 

-02,4.7726869E-03)  (0.6567292,-0.7539364) 

-0.7539366)  (6. 9561 149E-03, 1 . 5439859E-02) 

-02,4.0104548E-03)  (0.6131005,-0.7898440) 

-0.7898430)  <7. 7171684E-03, 1 . 3972549E-02) 

-02,3.3296335E-03)  (0.5675922,-0.8231732) 

-0.8231760)  ( 8. 304865 1 E-03, 1 . 2467677E-02) 

-02,2.731 3598E-03)  (0.5203235,-0.8538549) 

-0.8538533)  (8. 7 1871 41 E-03,  1 . 094455  IE-02) 


-02,2.2182921E-03) 


(0.4714787,-0.8817816) 
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pH 


(0.4714800,-0.8817813)  (8. 9602BB9E-03,9. 4234049E-03) 

15.31200 

(1. 18B1074E-02,1.791380BE-03)  (0. 4212167,-0.906B807) 

(0.4212141,-0.9068813)  (9. 0332972E-03 , 7 . 9225609E-03) 

15.40800 

(1.0937320E-02, 1 . 4509220E-03)  (0.3696674,-0.9290989) 

(0.3696675,-0.9290984)  (8 . 9446194E-03 , 6 . 4596324E-03) 

15.50400 

(9.9910377E-03, 1. 1954672E-03)  (0.3170410,-0.9483583) 

(0.3170418,-0.9483584)  ( 8 . 70 130 12E-03 ,5. 0533540E-03) 

15.60000 

<9.0487348E-03,1.0245877E-03)  (0.2634658,-0.9646257) 

(0.2634663,-0.9646257)  ( 8 . 3133159E-03 , 3. 71751 33E-03) 

15.69600 

(8. 1195952E-03,9.3581091E-04)  (0.2091443,-0.9778506) 

(0.2091438,-0.9778512)  <7. 7920314E-03.2. 4674323E-03) 

15.79200 

<7.2100945E-03,9.2572690E-04)  (0. 1542246,-0.9880092) 

(0. 1542242,-0.9880095)  (7. 1 49 1 729E-03 , 1 . 3154999E-03) 

15.88800 

<6. 3271616E-03,9. 9132408E-04)  (9 . 8894700E-02 , -0. 9950770) 

(9. 8896667E-02, -0.9950770)  (6. 3984897E-03 , 2. 7323177E-04) 

15.98400 

(5. 4778936E-03,1. 1289094E-03)  (4. 331 4267E-02,-0. 9990461 ) 

(4. 33 18227E-02, -0.9990457!  (5. 554871 2E-03 , -6. 507266  IE-04 ) 

16.08000 

(4.6668719E-03,1.3330920E-03)  (- 1 . 2329667E-02 , -0. 9999 1 24 ) 

(-1.2333964E-02, -0.9999121)  (4. 6325B24E-03, -1 . 4477773E-03) 

16. 17600 

(3. 9003449E-03, 1.5990 353E-03)  (-6. 7883283E-02 , -0. 9976846) 

(-6. 7881212E-02, -0.9976840)  (3. 6476655E-03 , -2. 1 1261 33E-03) 

16.27200 

(3. 1829933E-03,1.92153B0E-03)  (-0. 1231601,-0.9923795) 

(-0.1231597,-0.9923797)  (2. 6165834E-03 , -2. 641 4776E-03) 

16.36800 

<2. 51B5787E-03,2.2949139E— 03)  (-0. 1 78002B , -0. 9840239) 

(-0. 1780028,-0.9840238)  ( 1 . 5548397E-03 , -3. 031 B049E-03 ) 

16.46400 

(1.9108233E-03,2.7130966E-03)  (-0.2322484,-0.9726511) 

(-0.2322487,-0.9726509)  (4. 7B68307E-04 , -3. 2B37493E-03) 

16.56000 

( 1. 36261 1 4E-03, 3. 16961 18E-03)  (-0.2857049,-0.9583113) 

(-0.2857063,-0.9583109)  ( -5. 9565267E-04 , -3. 3983923E-03) 

16.65600 

(8.7576511E-04l3.65B'»052E-03)  (-0.3382418,-0.9410517) 

(-0.3382428,-0.9410510)  (-1 . 6540725E-03 , -3. 379 1540E-03) 

16.75200 

<4. 5196188E-04 ,4. 1737645E-03)  (-0.3896937,-0.9209350) 

(-0.3896935,-0.9209347)  (-2. 6812442E-03,-3. 23042 13E-03) 

16.84800 

(9. 2306254 E- 05, 4.7086221E-03)  (-0. 4399188,-0. B9B0247) 

(-0.4399184,-0.8980258)  (-3. 66385B6E-03 , -2. 9590100E-03) 

16.94400 

(-2 . 0387334E-04 ,5. 2566156E-03)  (-0.4887539,-0.8724062) 


m 
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(-0.4887477,-0.8724095)  ( 

17.04000 

(-4. 354 10658-04,5.81 182078-03) 
(-0.5360470,-0.8441678)  <■ 

17.  13600 

(-6.03997508-04,6. 3684662E-03) 
(-0.5816741,-0.8133960)  (■ 

17.23200 

(-7. 1176427E-04, 6. 91995168-03) 
(-0.6255079,-0.7801870)  (- 

17.32800 

(-7. 59441878-04,7. 4621327E-03) 
(-0.6674093,-0.7446533)  (- 

17.42400 

(-7. 50030278-04, 7. 9882778E-03) 
(-0.7072658,-0.7069025)  <- 

17.52000 

(-6.85818038-04,8. 49468 17E-03) 
(-0.7449560,-0.6670592)  <- 

17.61600 

(-5. 6950748E-04, 8. 97632458-03) 
(-0.7803739,-0.6252488)  (- 

17.71200 

(-4. 0656369E-04, 9. 4291 1028-03) 
(-0.8134203,-0.5815997)  (- 

17.80800 

(-1.9931207E-04.9.8498920E-03) 
(-0.8440095,-0.5362371)  (- 

17.90400 

(4.7SB191BE— 05, 1. 0234 153E-02) 
(-0.8720468,-0.4893157)  (- 

18.00000 

(3.30331B8E-04, 1 . 0580279E-02) 
(-0.8974622,-0.4409644)  (- 


-4. 589328 18-03, -2. 5712170E-03) 

(-0.5360438,-0.8441700) 
-5. 4451725E-03 , -2. 0775753E-03) 

(-0.5816740,-0.8133972) 
—6. 22167908-03,-1 . 4B7 25528-03) 

(-0.6255097,-0.7801852) 
—6. 9094244E-03,-8. 0975B61E-04) 

(-0.6674084,-0.7446541) 
-7.50031468-03,-5.91 45659E-05) 

(-0.7072648,-0.7069036) 
-7. 9878466E-03, 7. 54 104358-04) 

(-0.7449536,-0.6670628) 
-B.3676111E-03, 1 . 6159732E-03) 

(-0.7803729,-0.6252490) 

-8.6360276E-03,2.5132750E-03> 

(-0.8134216,-0.5815986) 

-8.7908590E-03,3.434346BE-03) 

(-0.8440076,-0.5362407) 

-B.8321501E-03,4.3649161E-03) 

(-0.8720467,-0.4893166) 

■8. 75964958-03,5. 2922373E-03) 

(-0.8974621,-0.4409642) 

■8. 5770553E-03 , 6. 2036063E-03) 
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APPENDIX  D 

COUPLING  COEFFICIENT  FOR  A  TEjq-45  DEGREE  TAPERED  RECTANCULAR  WAVEGUIDE 

The  following  derivation  uses  analytical  techniques  to  determine 
the  eigenvalue,  eigenfunction  and  coupling  coefficient  of  a  rectangular 
waveguide  having  a  constant  width  a,  and  tapered  height  b. 


Fig,  Dl.  TEjq  mode-45  degree  tapered  rectangular 
waveguide  with  constant  width. 


The  general  solution  to  the  Helmholtz  wave  equation 


$  ♦  ,  ,  +  h .  .  i|»f  .  *  0 

t  rtp]  [p]  ( p] 


for  TE  modes  having  the  boundary  condition 


(D.l) 
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$ 


0  on  C(x,y,z) 


(D.2) 


mxx  nxx 

ip.  ip.  ,  =  A  cos  — —  cos  —j— 

[pi  [tun]  mn  a  b 


(D.3) 


which  can  be  written  for  TF.jq  modes  as 


ho  ■  Aio cos  r 


(D.4) 


Eq .  D.4  is  the  eigenfunction,  and  its  eigenvalue  is  given  by 


10  a 


(D.5) 


The  unknown  normalization  constant  Ajq  can  be  determined  using  Solymar's 
normalization5 


where 


/  J  I 


e[P]l  da  "  1 


"(Pi  “  C[ 10]  =  Zo  X  V[10] 


(D.6) 


(D.  7 ) 


Differentiating  <J). 


V(10J 


.  X  ,  XX 
-A, .  —  sin  —  y 
10  a  a 


(D.8) 


*{8 


-  120  - 


and  substituting  the  result  into  Eq.  D.7  yields 


»  n  nx 

'no]  ■  -Aio  i sin  r y 


which  can  he  used  in  Eq.  D.6  to  solve  for  Ain  as 


(D.9) 


b  a  „  2  0 

/  /  A10  T  Sln  “  dx  dy 

0  0  1  a  3 


2  2 

b(z)ir  A  n  a  „ 

10  t  .  2  nx  . 

- = -  I  sin  —  dx 

2  a 

a  o 

2  2 

b(z)ir  A1q 


(D.10) 


l  (  2a  V 


(D.ll) 


According  to  Appendix  A,  the  backward  coupling  coefficient  for  the  TE, 


mode  is 


’[10][10] 


-U 


tan  9 


D.12) 


C(x,y,z) 


-  < 


Since  only  the  top  and  bottom  of  the  waveguide  are  tapered,  the  sides  do 
not  contribute  to  the  contour  integral  and  Eq.  D.12  can  be  rewritten  as 
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'[10] [10] 


-  -  —  /  tan  0 


m 


al  2 

r  1  2a  ^  0  ir  .  2  ttx  , 

1  t  fur  6  t  sin  r dx 

(j  it  a 


tan  6 
b(z) 


(D.13) 


The  rule  for  the  sign  of  tan  0  is:  walls  that  taper  towards  the 
z-axis  have  negative  0  ,  those  that  taper  away  have  positive  0.  This 
stems  from  the  following  expression  for  tan  0, 


tan  0 


b2  -  bl 


(D.14) 


Substituting  Eq.  D.14  into  Eq.  D.13  gives  the  final  expression  for  the 
TE^q  backward  coupling  coefficient  of  a  rectangular  waveguide  having 
constant  width  and  tapered  from  bl  to  b2. 


bl  -  b2 

[ 10 ] [10]  “  2  b(z)  L 


(D.15) 


In  Section  III,  Eq .  D.15  is  used  as  the  standard  against  which  computed 
results  were  compared.  The  test  case  is  that  of  a  waveguide  like  the 
one  shown  in  Fig.  D.l,  having  a  taper  of  0  *  -45°.  The  coupling  coeffi¬ 
cient  can  be  directly  computed  from  Eq .  D.15  as 


’[10] [10]  "  b(z) 


(D. 16) 


mmmm 


A  general  equation  for  the  linear  variation  of  waveguide  height  with 
axial  position  z  is 


b(z )  -  bl  ~-'Z)  +  b2  £ 


(D.17) 


In  the  example  of  Section  III,  the  cross  section  referred  to  could  be  at 
any  position  along  the  z-axis.  For  the  case  of  this  example,  z  =  0  and 


b(z)  *  bl 


(D. 18) 


'[10] [10]  *  bl 


(D.19) 


The  results  of  changing  the  mesh  size  of  a  (10  x  20)  cm  rectangular 
waveguide  are  shown  in  Table  1.  Since  bl  =  10.0  for  this  case, 


S[10] [10]  "  0,1  cm 


(D.20) 


It  is  clear  that  the  computer  program  accurately  calculates  the  coupling 
coefficient  of  nonridged  waveguides  to  within  ~  0.3  percent. 
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APPENDIX  E 


1 

’flJj 


FIELD  NORMALIZATION 


In  order  to  obtain  a  correct  value  for  the  coupling  coefficient 


S [ i © ] [io]’  tlie  numerically  obtained  eigenfunction  must  be  properly 


normalized.  The  proper  normalization  equation  is  given  in  Section  II  by 


Eq .  6 . 


A(x,y,z)  P 


*  2 

e  da  =  1 


(E.l) 


By  using  standard  vector  identities,  the  two-dimensional  form  of  Green's 


Theorem1 7 


/  - 
A 


+2  >> 

<P  <p)  da 


'c 


♦  lx)  da 


(E.2) 


the  Helmholtz  Wave  Equation 


+2  2 
V  <p  +  h  ip  =  0 

t  p  p  p 


and  either  of  Eqs.  4  or  5  with  their  appropriate  boundary  conditions,  it 


can  be  shown  that  Eq .  E.l  reduces  to 


I  l*p|2  da  -ij 

A  h 

P 


(E.3) 


The  numerically  obtained  field  configuration  ip  can  be  normalized  by 

pN 


integrating  its  magnitude  squared  over  the  guide  cross  section  A  with 
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VI  -.-'V  OVVV'-.-VV'. 


t*ht 


the  result 


/a  l*pNl  <*“  -  » 


In  order  to  match  the  field  normalizations,  it  is  necessary  that 


(E.4) 


* 


4*. 


pN 


p 


1/2 


(E.5) 


The  finite-difference  algorithm  in  the  ridged  waveguide  program  con¬ 
verges  to  the  value 


2  2 

u  h  =  X 
P  P 


(E.6) 


thus 


(X  ) 
_P_ 


1/2 


(E.7) 


Since  up  is  the  numerical  approximation  of  Solymar's  hp,  Eq .  E.7  can  be 
substituted  into  Eq .  E.5  and 


=  t 


pN 


(»‘P) 


1/2 


(E.8) 


Using  Eq .  E.8,  every  numerically  obtained  field  point  can  be  scaled  to 
Solymar's  normalization.  This  puts  the  mode  eigenfunction  in  the  appro¬ 
priate  form  for  computing  Solymar's  coupling  coefficients. 
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